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THE RISK OF SOLAR PROTON 
EVENTS TO SPACE MISSIONS 

SUMMARY 


The total dose from solar protons was tabulated in rads-tissue for 
weekly time intervals, and the number of weeks that gave a dose above 25 rads 1 
behind 10 g/cm 2 of aluminum for the active six years of the 19th cycle were 
called dangerous or large event weeks. Only three such event weeks were 
found during the past 20 years. 

Even though the possibility of smaller events was examined, it was 
found that for any reasonably high confidence level ( 95 percent) , the smaller 
events could be ignored. Consequently , the total particle flux for the 1 9th 
cycle was divided by a factor of three , and a single large event week was 
determined. 

Using this spectrum, one can calculate the tissue dose in rads at the 
center of an aluminum spherical shell. To correct for geometric effects and 
self-shielding, this dose should be reduced by a factor of approximately three. 

The Poisson distribution seemed to be the most logical choice for pre- 
dicting the probability of an event occurring. Section III is devoted to examin- 
ing this conclusion and the methods of using this probability model. Section 
IV defines the confidence one could use in employing the Poisson process and 
specifically arriving at confidence levels for the experimental or observed 
value of the mean of the Poisson distribution function. Several examples are 
given for different mission lengths, and comparison to the results of other 
authors are made. 

Section V is an extension of the Poisson process that incorporates the 
concept of small sample theory and arrives at the expected distribution function 
that answers the following question. If events are observed in time T^, 

what is the probability of seeing X events in any observation time t? 


i. 


To convert rads to SI Units in j oules/kilogram , multiply rads by 

0 . 01 . 



SECTION I. INTRODUCTION 


This study is concerned with the practical treatment of hazards result- 
ing from solar proton events outside the magnetosphere of the earth. It is not 
concerned with prediction of flares as such or the long-range solar cycle in- 
dicators such as sunspots. This distinction is made because these phenomena 
have been related to solar proton events and it has been assumed that statistical 
studies of these phenomena give insight into the problem of dangerous solar 
proton events. In fact, the study of all solar proton events as such is only 
casually related to the problem of large dose rates inside realistic spacecraft. 
For example, in Reference 1, 76 events are listed from 1942 through 1963. 
Since 1963 there are probably an additional 24, which yields a statistical 
sample of some 100 events. However, the correlation of this large sample 
with the events that are of real danger to space flight is very poor. For ex- 
ample, only six large events from 1950 through 1969 would have given about 
85 percent of the total 20-year proton dose behind a thin wall of only 2 cm 
of aluminum. In addition, these six flares occurred in three weekly periods. 
Thus, three of these six large events occurred from July 10 through July 16, 
1959; two occurred from November 12 through November 18, 1960; and on 
February 23, 1956, one large event was observed. The duration of an event 
is from 1 to 3 days. Even with this small sample it seems that if conditions 
are suitable for a large event, the chances are very good that it will be followed 
within hours by another large event. Thus, the treating of such large events 
as mutually independent seems questionable. (For example, if three tornadoes 
hit a city within a 24-hour time period and this occurred once in 30 years , one 
would hesitate to say that the chance for a tornado to hit this city is one in 10 
years. ) Perhaps a time span of at least 1 week should be used to depict a 
total solar event; it would be designated as the solar proton flux or dose per 
week. Under this definition, there were only 3 sample weeks of large solar 
events from about 1950 through 1969 ( 1000 weeks) . One is not dealing with an 
ordinary problem of statistical analysis, but with rare large events. 

The author has undertaken the awesome task of predicting the improb- 
able not by virtue of his background in solar physics or statistical analysis 
but because of his concern with protection of man and his radiation-sensitive 
equipment from space radiations in applications involving realistic space- 
craft and missions. However, the reader should ask, "Why not leave this 
field of statistical astronomy to the experts?" The answer is that the environ- 
mental scientists do not have to design or evaluate realistic shields for sensi- 
tive film or radiation-conscious astronauts. 
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The consequence of having several solar proton prediction models 
(which the author does not wish to evaluate) has led to a wide disparity in re- 
sults , especially when a reasonably high level of statistical certainty is desired. 
For example, at the 99-percent probability (percentile) level, the predicted 
dose behind 20 g/cm 2 of aluminum for a i-year mission may vary by a factor of 
10 or more between different investigators [ 2] . This may not sound too bad 
considering the nature of the problem; however, if a mission is planned with the 
requirement that the astronauts should not receive doses exceeding a 100-rad 
skin dose with a 99-percent probability and if one solar proton dose prediction 
model requires a shield of only 15 g/cm 2 while another model requires a 
shield of 50 g/cm 2 , the radiation analyst has to make a vital decision that 
possibly affects the life of an astronaut. Since no one wishes to be responsible 
for making a decision that could lead to dire consequences, the most pessi- 
mistic model is often chosen. The author does not question this approach ex- 
cept to note that the desire to be on the safe side may readily become uncon- 
trollable and a subsequent loss in mission capability may develop. 

As a result of the problems mentioned above, NASA has been an easy 
target for various academic and industrial contractors who propose to solve 
certain aspects of this problem. Their proposals may range from early solar 
flare predictions based on observing the sun to statistical studies of sunspots 
and recorded flare data , and include special concepts for shielding against 
large solar proton events. Many of these studies are interesting and even of 
credible workmanship, but a large majority are of only marginal value to a 
radiation analyst who is confronted with real problems. This writer feels 
that most of the studies are of dubious value , because statistical studies of 
sunspots or flares or all solar proton events do not provide any meaningful 
guide to predicting large, very hazardous proton events. For example, there 
were only 3 weeks from the past 1000 weeks during which an astronaut inside 
a spacecraft having a wall thickness of only 3 cm of aluminum would have re- 
ceived as much as a 20- to 30-rad dose at a depth of 5 cm in the body. 

In addition to the statistics problem, there exists the rather common 
procedure of calculating the radiation hazard from a solar proton event for a 
point detector at the center of a spherical shell. This may lead to overesti- 
mates of the actual dose to an astronaut in a realistic spacecraft by a factor 
of from approximately three to as much as six. Also, it should be mentioned 
that the solar proton dose inside the earth's magnetic field will be consider- 
ably reduced. The actual reduction depends on the orbital inclination and 
altitude of the spacecraft, as well as the disturbances in the earth's magnetic 
field that are caused by magnetic storms accompanying the solar proton event. 
Treating this aspect of the problem correctly requires a sophisticated mathe- 
matical model; for an example, see Reference 3. 
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SECTION II. A PRIMITIVE DOSE MODEL IS PROPOSED 


The most direct measure of the hazard of a given solar proton event is 
the rads-tissue absorbed dose that would be measured behind various thick- 
nesses of a typical spacecraft material such as aluminum. The simplest 
method is to find the point tissue dose at the center of a spherical shell of 
aluminum. However, because of self-shielding by the astronaut (approximate- 
ly a factor of two, the spacecraft geometry, and onboard equipment, the actual 
skin dose may be less , by a factor of three or more , than the point dose at the 
center of a spherical shell. This factor varies depending on the solar proton 
spectrum, the spacecraft geometry, and the location of the astronaut in the 
craft. This writer suggests using a reduction factor of approximately three 
for the solar proton point dose at the center of a spherical shell to estimate 
the probable skin dose to an astronaut in a real spacecraft of a given average 
thickness. This factor is believed to be reasonable in the sense that this 
reduced dose is probably high when the complex geometry and spectrum shape 
of a solar proton event are considered. 2 To correct for depth dose (bone- 
marrow depth) , a thickness of 5 cm of tissue is often employed. This is 
approximately equal to 6. 5 g/cm 2 of aluminum in equivalent shielding effective- 
ness. 


To clarify the exact assumptions that are to be employed, the following 
information is pertinent. The only adequate data available at this writing for 
solar proton predictions are from the 19th cycle (1954 through 1964). However, 
there have been observations of sunspots (indicators of solar activity) for 
approximately 200 years. Based on more than 200 years of observations, in 
which the first cycle dates back to the middle of the 18th century (the average 
cycle length is about 11 years) , the 19th cycle had the highest maximum 
sunspot count yet recorded. In fact, the average maximum sunspot count is 
more than a factor of two lower than that of the 19th cycle. The 20th cycle, 
which we are now well into (past the peak activity) , has a sunspot count 
somewhat above this average. There has not been a large solar proton event 
in the 20th cycle comparable to the eight largest events of the 19th cycle 
(dose behind 10 g/cm 2 ) . At this point one may be led to believe that the 19th 
cycle is a relatively rare cycle. With the present low occurrence of large 
solar proton events and the rather extensive sunspot counting dating back over 
200 years, one might conclude that the probability of a solar cycle being as 
active as the 19th is on the order of 1/20 or 0. 05. This, of course, cannot be 


2. It should be noted that this factor of three may be too large for a very 
"hard” spectrum such as seen in the trapped radiation belts , where a 
a factor of 2. 5 may be more appropriate. 
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objectively demonstrated, and will not be, unless considerably more knowledge 
is obtained about the physics of the sun. One valid objection to the above is 
that the number of sunspots is a poor indicator of large solar proton events. 
Also, the sunspot indicator may have changed during the last 50 years because 
of better observation techniques, etc., so that possibly the first 15 of the 20 
observed cycles should not be used. The purpose of this paper is not to attempt 
to evaluate the above but to present the information for the reader's considera- 
tion. At any rate , the reader may not find it difficult to believe that the use 
of the 19th cycle solar proton flux data may yield a pessimistic estimate of the 
solar proton hazard. With the above background, one can at least study the 
proton events characterized by the 19th cycle and infer proton events for future 
cycles that are similar to the 19th cycle. 

Next, consider the observation that during the 19th cycle, the large 
solar proton events did occur primarily around the most active years of the 
11 -year cycle, however not necessarily in proportion to sunspot count. It is 
generally assumed that there were only about 6 years of observed large solar 
proton events. Thus, the following analysis is based on the so-called active 
6 years (300 weeks) of the cycle. The second assumption is that it is fair 
and logical to lump the actual dose rates over an active week into units of total 
dose per week (see Section I). This may include up to three proton events 
in a given week. Thus , instead of using fundamental units of proton events , 
it is proposed that observed weekly dose rates behind various aluminum thick- 
nesses be used. The data available from various sources gave a total of 24 
proton events that were worth considering. ( The ten or more weeks discarded 
gave less than 3 rads total behind 5 g/cm 2 of aluminum. ) The results of 
grouping this data into weekly time periods gave a total of 18 weeks with the 
following frequencies: 13 weeks, one event; 4 weeks, two events; and 1 week, 
three events. Most of this data is recorded in NASA TN D-4404 [ 3] . The 
only major revision in Reference 3 is in the use of the spectrum of A. J. 

Masley [4] for the November 12, i960, flare. The previous calculations 
were based on the work of W. R. Webber [ 5] . The differences are shown in 
Table 1. 

One note of explanation should be made regarding the dose rates that 
are used in this study. They include a correction for secondary particles 
and thus are the sum of the primary proton dose and the secondary particles 
(neutrons and protons) . This total dose at 20 g/cm 2 and at 10 g/cm 2 is 
approximately 20 percent and approximately 10 percent, respectively, above 
the primary proton dose rates. 

Now one comes to a very treacherous part of the analysis. How does 
one choose a rational and sufficient model for an event dose-week for solar 
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TABLE 1. SKIN DOSES (RADS- TISSUE) BEHIND VARIOUS SHIELDS 
OF ALUMINUM FOR NOVEMBER 12, 1960, FLARE 


Shield 



^X^tg/cm 2 ) 



Source 

2 

5 

Webber 

270 

106 

Masley 

1900 

400 

i_ 


a. 20 g/cxn 2 -aluminum = 7. 4 cm = 2. 9 in. 


10 

45 

61 


a 

20 

16. 2 
9.5 


protons ? Table 2 is presented as a summary of the weekly solar proton doses 
of the 19th solar cycle that the writer will consider. The absorbed doses in 
rads-tissue were calculated at the center point of a spherical shield of aluminum, 
with the shell thickness being designated as the shield. 

The data of Table 2 have been grouped into three categories: (1) large 
weekly doses, (2) medium weekly doses , and (3) small weekly doses. Cate- 
gory 1 , consisting of only 3 weeks , gave more than 84 percent of the total 
6-year dose. Category 2 with 6 samples gave about 13 percent, and category 
3 with 9 samples gave less than 3 percent of the total dose. 

The grouping shown is certainly not unique , and the last 2 weeks of 
category 2 perhaps should be in category 3, or perhaps the largest dose week 
of category 2 should be in category 1. The following work has attempted to 
provide for various combinations that the reader may wish to investigate by 
making the methods of approach more important than the choice of a precise 
set of data. 

Since over 84 percent of the total solar proton dose is grouped into the 
3 weeks of the first category, it would seem that one feasible solar proton 
dose model for an active week could be depicted by dividing the total 6-year 
dose behind the various shields by three. Since the most active phase of the 
19th cycle is about 6 years, it seems sufficient to use 300 weeks as the basic 
period of major solar activity. With this assumption, which will be used 
throughout this study, it follows that the estimated weekly expectation of a 
dangerous large solar proton dose would be = 0. 01 ( event/ week) . Thus, 
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TABLE 2. WEEKLY TISSUE DOSES OF THE 19TH SOLAR CYCLE 


Sample 



Number 

of 

Shield 

(6-Year 

Dose) 


Number 

Week 

Year 

Flares 

5 ( 906) a 

10(251) 

20(73) 

Remarks 

1 

Nov 12 - 18 

1960 

2 

456 

83. 4 

17. 0 

Category 1 — Large Weekly Doses 

2 

July 10 — 16 

1959 

3 

216 

75. 0 

21. 2 

X ^ = 0. 01 

3 

Feb. 23 

Summed Dose 
Percent of To- 
tal 6 -Year Dose 

1956 

1 

(6) 

92 

764 
84. 3 

50. 0 

208. 4 
83 

24.8 

63. 0 
86. 3 

Approximately 85 Percent of 6-Year Total 

1 

May 10 

1959 

1 

59. 3 

18. 3 

4. 4 

Category 2 - Medium Weekly Doses 
X * 0. 02 

2 

July 18-20 

1961 

2 

22. 0 

8. 2 

2. 4 

2 

3 

March 23 

1958 

1 

10. 9 

2. 5 

0. 4 


4 

July 7 

1958 

1 

10. 5 

2. 3 

0. 4 


5 

Jan. 20 

1957 

1 

8. 3 

1. 8 

0. 3 


6 

Oct. 20 

Summed Dose 
Percent of To- 
tal 6 -Year Dose 

1957 

1 

(7) 

4. 1 

115. 1 
12. 7 

1. 8 

34. 9 
13. 9 

0. 7 

8. 6 

11. 8 

Approximately 13 Percent of 6-Year Total 

1 

Aug. 22 — 26 

1958 

2 

5. 9 

1. 0 

0. 2 

Category 3 — Small Weekly Doses 

2 

Aug. 29 

1957 

1 

4. 2 

0. 8 

0. 1 

X 3 " 0. 03 

3 

Nov. 20 

1960 

1 

3. 6 

1. 5 

0. 05 


4 

Aug. 3 

1956 

1 

2. 2 

1. 0 

0. 4 


5 

Sept. 13 

1960 

1 

2. 9 

1. 2 

0. 5 


6 

Aug. 1G 

1958 

1 

1. 8 

0. 4 

0. 1 


7 

July 3 

1957 

1 

1. 2 

0. 43 

0. 03 


s 

July 12 

1961 

1 

1. 4 

0. 30 

0. 04 


9 

April 28 — 29 

Summed Dose 
Percent of To- 
tal 6-Year Dose 

1960 

2 

(H) 

0. 77 

23. 97 
2. 65 

0. 27 

6. 90 
2. 75 

0 

1. 42 
1. 95 

Approximately 2 Percent of 6-Year Total 


a. The first number is the shield thickness in g/cm 2 -aluminum; 

the number in parenthesis is the 6-year dose behind the shield. 
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I III IHIII III llll 


the simple expectation for 300 weeks is three events, where each event yields 
one-third of the 19th solar cycle dose behind various shields. This concept 
is usable, but as the reader might suspect, the calculation of any fairly sophis- 
ticated statistical measure would lead to very large uncertainties because of 
the uncertainty in the value of the parameter (A Q ), since the actual number of 

recorded events is very small. Also, it might seem that this concept could 
not be used for predicting the likelihood of solar proton doses for short ex- 
posure periods. In the next section this point will be discussed in more detail. 
In any case , for a short exposure time , there is a greater chance for an event 
of category 2 or category 3 to occur than there is for an event of category 1 
to occur. For this reason a model will also be developed to reflect medium 
and small flare doses. Thus, two solar proton dose models will be constructed, 
and the results will be compared in the following work. 

The next part of our approach is to choose the best composite 19th 
cycle solar proton flare data that are available and to arrive at total tissue 
dose in rads from primary and secondary particles at the center of spherical 
shells with varying thicknesses. T. T. White, et al. [6] have developed a 
composite model (designated as the MSC model) for the total proton flux dur- 
ing the 19th cycle. This model gives a larger dose for energies less than 
approximately 115 MeV 3 * than does a model based on Webber's work [5]. 

Above 115 MeV a composite flare model based on Webber's work (designated 
as the MSFC model) gives a larger dose. The 19th-cycle, composite, 6-year 
proton spectrum that will be used in this work is: 

Model I (MSC) J(>p) = 5. 28 x 10 11 e~ P//73 (30 £ E < 115 MeV) 

, ( 1 ) 

Model II (MSFC) J(>p) = 1. 14 x 10 11 e _P//100 (E > 115 MeV) 

where p is in rigidity units of megavolts and J(>p) is the integral spectrum 
(proton/cm 2 ) with energies above p . For prot ons, the relationship between 
p(MV) and E(MeV) is given by p = V E 2 + 1876E [3] . 

Using the spectrum of equation (1) , the best estimate of the 6-year 
(300-week) total dose (primaries and secondaries) is given in Figure 1. From 
Figure 1 and Table 2 the solar proton 1-week-dose event models are con- 
structed. (Figure 1 gives the magnitude and Table 2 gives the fractions. ) 

Model I will be the dangerous solar event model that is represented by three 
large events over a period of 300 weeks. The expected number of events 


3. To convert electron volts to SI Units in joules, multiply electron 

volts by 1. 60210 x 10~ 19 . 
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RADS - TISSUE 



Figure i. Total tissue dose at the center of a spherical shield from solar 
protons of the 19th cycle as a function of shield thickness X. 
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in 1 week is given by A^ = 0. 01 (event/week) . Model II will depict the possi- 
bility of a large, medium, or small dose in a week where the percentages of 
the total dose in Table 2 are used to determine the relative size of a dose week 
in the three groups. For the three categories of Table 2, the values A^ = 0.01, 

A = 0. 02, and A = 0. 03 are the mutually independent weekly expectations of 

A o 

large, medium, and small doses, respectively. Table 3 summarizes the dose 
models. It should be clear that the author has presented only the primitive 
elements of a probability model in Table 3. Thus , the table gives the small 
sample estimates of the weekly expectations (A), and the consequences (rads/ 
week) if an event occurs according to model I or model II. It should also be 
made clear that the values of A are valid only for the 300 most active weeks 
of a solar cycle similar to the 19th. The major unknown factor is an estimate 
of the chance of obtaining a solar cycle that would give total proton doses as 
large or larger than the 19th cycle. 

The values of A should be considerably smaller for the so-called quiet 
sun (260 weeks) or else the magnitude of an event should be reduced to the 
level of an event of category 3. Perhaps one might assume that A = 0. 02 
(one event per year) during the 260 weeks of the quiet sun, but an event should 
be depicted by the D column (small event) of model II. A similar set of 

o 

revisions might be in order if the solar activity is definitely smaller than the 
19th cycle at the maximum point of the cycle. However, there is no justifica- 
tion for any of the above suggestions. 

As a conclusion to this section, it is of interest to recall that Figure 1 
represents the total dose versus spherical shell thickness for the 19th cycle. 

If for the data of Figure 1, an effort is made to correct for the self-shielding 
of an astronaut and the complex geometry of a spacecraft that has aluminum 
walls of 13. 5 g/cm 2 (assume that a 5-cm depth dose corresponds to an addi- 
tional aluminum thickness of 6. 5 g/cm 2 ) , then the 5-cm depth dose is esti- 
mated by dividing the dose at 20 g/cm 2 by a factor of three. Thus, one finds 
that the total dose from the 19th cycle was 25 rads at a 5-cm tissue depth. 

Of course, this value is for the solar proton dose and does not account for the 
galactic cosmic ray dose that may range from 5 to 12 rads/year behind 20 
g/cm 2 , depending on how much dose is contributed by the heavy cosmic ray 
particles ( Z > 2) . The implications of all the above are simply that since the 
19th solar cycle was an example of a very active sun and if one believes that 
the high energy proton fluxes received during this cycle are associated with 
this activity, he must be cautious in drawing certain conclusions. For ex- 
ample, the same astronaut in the previous example may have received from 
30 to 70 rads from galactic cosmic rays during a six -year trip, whereas the 
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TABLE 3. MODELS OF A i-WEEK DOSE EVENT FOR THE ACTIVE 
300 WEEKS OF THE 19TH SOLAR CYCLE 


Model I 



Model II 




N = 3 a 
100 Percent 
A Q = 0. 01° 


N( Large) = 3 
84 Percent 
= 0. 01 

N ( Med. ) = 6 
13 Percent 
A g = 0. 02 

N( Small) = 9 
3 Percent 
A = 0. 03 

d 

Z D i N i 

i=l 

X 

(g/dm 2 ) 

/ rads \ 
oVweek/ 

X 

(g/cm 2 ) 

D ( rac * s \ 
i\week/ 

D ,( rad ?) 
2 Vweek/ 

D,( ra t) 

3 Vweek/ 

Total 

Dose 

300 

Weeks 

2 

1420 

2 

1195 

95 

15 

4290 

5 

340 

5 

286 

24 

3.4 

1032.6 

10 

93 

10 

78 

7. 5 

0. 9 

281. 1 

15 

43 

15 

36 

3. 5 

0. 43 

132. 87 

20 

25 

20 

21 

2. 0 

0. 25 

77. 25 

30 

11 

30 

9. 3 

0. 80 

0. 11 

33. 69 

40 

6 

40 

5. 03 

0. 45 

0. 06 

18. 33 

60 

2 

60 

1. 70 

0. 15 

0. 02 

6. 18 


a. N indicates the number of weeks used in the dose category. 

b. Percentage shown is the percent of total dose in the category. 

c. A gives the expectations (probability) of the category for 1, week. 

previously discussed dose from solar protons at bone-marrow depth was 25 
rads (13. 5 g/cm 2 aluminum walls) . Here, one assumes that the galactic 
cosmic ray spectrum is so energetic that geometry factors are negligible for 
dose reduction. Thus, the cosmic ray dose component may determine the 
limiting dose factor for long duration space travel. The purpose of the fore- 
going discussion is not to minimize the importance of solar proton events , but 
to illustrate that it may be quite feasible to shield against solar protons but 
probably impractical to consider shielding against galactic cosmic rays. This 
very fact may be more important in determining man's exposure time to space 
outside the earth's magnetic field than the solar protons. 
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SECTION III. A PROBABILITY MODEL IS DERIVED 


In this section, the problem of choosing and using a probability density 
function in order to arrive at the probability of getting X events in t weeks 
and the associated problem of probability or percentile levels shall be ad- 
dressed. The important question of establishing confidence intervals for the 
basic statistical parameters that are obtained from a small sample will be 
undertaken in the next section. In this section, the naive assumption that the 
basic statistics or population parameters (mean and variance) are well known, 
either by experience or a priori knowledge, will be made. 

To derive the probability model that seems to be the most natural out- 
growth of the solar proton dose week, the process will be described in terms 
of probabilities P^(t) that exactly n events occur during a time interval t 

(weeks in our case). Thus, P^(t) is the probability of no event in the interval 

t, and 1 - P^(t) is the probability of one or more events. Next, A is defined 

to be the mean or expectation of an event for a unit time interval. That is, 


A 

X 


number of ev ents 
total weeks of observation 


( 2 ) 


is a statistical estimate of X. More precisely, X is a constant that deter- 
mines the density of points on the t axis. Thus, for a small interval of time 
At, the probability of one or more events is given by 


i - P (At) = A At + e (At) 


(3) 


where e (At) is an infinitesimal and small compared to A At such that 
€ (At) 


lim 

At-*0 


At 


= 0 


Now, the following postulate is made: 

Whatever the number of events during (0,t) , the probability during 
(t, t + At) that one event occurs is given by P^(At) = A At + (At) , 

and the probability that more than one event occurs is given by 

P n>l (AW ' £ l< A «- 
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These conditions are the basic assumptions of the Poisson process 
[7] . It should be clear that for a small time interval At, the chance for one 
event is approximately AAt and the chance for more than one event is very 
small compared to A. At. Since t is for a relative time scale, it is not con- 
tradictory that a At of 1 week can be small on our time scale. The above 
conditions lead to a system of differential equations for P q ( t) . Consider two 

time intervals (0,t) and (t, t + At) , where At is small . If n ^ i, n events 
can occur in the interval (0, t + At) in three ways: 

1. No events during (t, t + At) and n events during (0, t). 

2. One event durihg (t, t + At) and n-i events during (0, t). 

3. X 2: 2 events during (t, t + At) and n-X events during (0, t). 

Following our assumptions , the probability for the first case is given by 

P (t) P (At) = P (t) (1 - A At) + e . (At) , (4) 

n 0 n l 

where P^(At) = 1 - A At - (At) from equation (3). 


For the second case, one has similarly 

P ,(t) P, (At) = P (t)AAt+ € (At) 
n- 1 l n-i A 


and for the third case , 


P X = 2 (t > p n-x (it > = V At) 


The result of combining the above is 

P (t + At) = P (t) (1 - A At) + P At) A At + e'(At) 

n n n-1 


which can be expanded and rewritten as 


P (t + At) - P (t) , 

“ — ^ — = - AP (t) + AP (t) + 

. , n n-1 At 

At 


(5) 


( 6 ) 


(7) 


( 8 ) 
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As At-*-0 the last term tends to zero, hence one obtains in the limit 


a? «) 

— -2— = - XP (t) + XP ,<t), na 1 . (9) 

at n n— 1 

For n = 0, the second and third possibilities above do not exist, and equation 
(7) becomes simply 

P Q (t+ At) = P Q (t) (1 - A At) + e Q (At) , (10) 


which leads to 


P Q (t + At) - P Q (t) 
At 


- XP 0 (t > 


e o (At) 

At 


(ID 


and in the limit 


diyt) 

dt 


- xp o (t) 


( 12 ) 


From equation (12) and P^(0) = 1, one obtains 


p 0 (t> 


-At 


(13) 


Using equation (13) and P^(0) = 0, equation (9) can be solved for P^(t) = Ate 
Using the fact that P^O) = 0 (n>0) » equation (9) becomes a recursion equation 
and successive values of P (t) can be found. The resulting solutions give the 
terms of the Poisson distribution: 

n 


-At 


P (t) = e" 
n 


At (At) 


n : 


(14) 


If the reader accepts the postulate following equation (3) , the Poisson 
distribution is the natural outcome. The foregoing arguments have been pre- 
sented to minimize the illusion that the author has "pulled a distribution 
function out of the sky. " For a more rigorous treatment, Reference 7 is 
recommended. However, there is one condition necessary for a Poisson pro- 
cess that can be examined easily. This condition is that the occurrence of 
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events on the time axis (0,T) are random, independent, and uniformly dis- 
tributed. That is, no time interval At in (0,T) is more likely to have an 
event than any other time interval of the same length. To test this condition 
for the data of Table 2, the following development is given. First, the time 
interval (0, T) of our solar proton events is considered to be 300 weeks. Next, 
the time from the starting point zero to the occurrence of the ith event is de- 
fined to be t. and the sum of n such events in (0,T) is defined to be 


S n - i \ . <«> 

ItI 

Now, for moderately large values of n, the sum of n independent random 

variables , each uniformly distributed on the interval 0 to T , may be considered 
to be normally distributed with mean, 

E ( S n) ' “(i) ’ < 16) 

and variance, 

Var ( s „)= “(tf) < 17 » 

The standard deviation is then 

°( S n)- 7t(f) (18) 

(See Reference 8, Stochastic Processes by E. Parzen. ) The reader may 
object to the use of this test since there is only one set of sums S^ to test, 

and the values of n are small for either large, medium, or small events 
(n = 3, 6, 9). In any case, the test should have some merit in establishing 
a degree of credibility. To carry out the above test, the time scale will be- 
gin on January 1, 1956, and terminate on October 7, 1961. This time span 
corresponds to 300 weeks. To aid in showing any tendency for lumping, a 
further division of the event occurrence time into three consecutive groups 
of 100 weeks is made. Table 4 summarizes the results; the order of events 
is identical to that of Table 2. The results of Table 4 tend to imply that the 
assumption of random, uniformly distributed events on the time axis is valid 
at a reasonably high level of certainty. The assumptions of a normal dis- 
tribution in for moderately large n could allow a probability statement 
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to be made concerning the degree of belief that the distribution is random and 

uniform; however, since n is rather small, it is believed that the value of 

such a statement is dubious. Note that all values of the calculated S fall well 

n 

within a one-sigma bound of the theoretical mean. 

There are many interesting uses of the Poisson distribution in addition 
to the occurrence of rare events in a continuum of time. The distribution is 
used to approximate the binomial distribution for the case of rare events 
(p < 0. 05) . The word rare means individually rare. In a large population 
several such events may occur, but the probability of occurrence of each in- 
dividual event is small; for example, the number of people killed by horses in 
1969. An important feature of the Poisson distribution is that for large values 
of the mean ( At » 20) , the distribution approaches the normal (or Gaussian) 
distribution. There are many applications of the Poisson distribution given 
in any standard text on probability and statistics. The most common include 
such studies as the number born blind in a large city, radioactive disintegra- 
tion, bacteria on plates, telephone traffic, etc. The remainder of this section 
will be devoted to a discussion of the Poisson distribution and how to apply it 
to our class of problems. 

If a discrete variable X has a Poisson distribution, then the probability 
that (X = x) is given by the expression: 

-m x 

P (X = x) = 6 v , m , x = 0, 1, 2, . . . ; (19) 

and P(X = 0) = e m , since 01 =1. The mean (or expected value of the 
variable X) is given by n = m, and the variance a 2 is given by cr 2 = m. The 
standard deviation is given simply by cr = \Tin. The summation equation, 

Yj e~ m m X /x ! = 1 , is satisfied. Since the time -dependent form of the 
x=0 

Poisson is primarily of interest, m = At becomes the mean. Even though the 
Poisson is discrete in the variable X, it is continuous in t. For purposes 
herein, the Poisson distribution shall be written as: 

-At , ,x 

P (X = x) = , x = 0, 1, 2 (20) 

X • 

where t = the number of weeks, x = the number of events in t weeks, A = the 
mean number of events per week with At becoming the expected or mean 


M439 


17 



number of events in t weeks, and P(X = x) is the probability of exactly x 
events in t weeks. 


Thus, 


P(X2t x*) = ^ 

x=x* 


-At,. .. x 
e (At) 


( 21 ) 


and 


x* -At ,x 
e (At) 


P(Xs x *) = 2 
x=0 

Note that since one is dealing with a discrete variable x, 

P(X >x*) = P(X ?= x*) - P(X = x*) = 1 - P(X =£ x*) 


( 22 ) 


(23) 


and 


P(X < x*) = P(X < x*) - P(X = x*) =l-P(X^x*) . (24) 

Tables of the Poisson distribution are available in most handbooks of statistics. 
They are usually given in the form of the discrete distribution of equation (19) 
and the summed Poisson P(X ^ x*) . Other combinations can be obtained by 
using the relations shown in equations (23) and (24). Figure 2 depicts a 
family of summed Poisson curves as a function of At, obtained by cross plot- 
ting from a standard table. For example, if A = 0. 02 and t = 100, then 
At = 2. From Figure 2, the probability of getting six or more events is approx- 
imately 0. 0166 and the probability of getting less than six (five or less ) events 
is 0. 9834 [from equation (24)]. 

If one wished to find a value of X that would not be exceeded, for 
example, at least 99 percent of the time (a i-percent chance at most of ex- 
ceeding X) for a given At, the cumulative distribution function P(X ^ x*) is 
utilized. The value of X that is sought here is the smallest x* that satisfies 


Pr (X =£ x * ) = 

x* -At^.-x 

Z * p 

f 

(25) 


x=0 



where P = 0. 99 for the 

case in question. 




18 


PROBABILITY THAT X > K EVENTS AS A FUNCTION OF At 


1.0 



Figure 2. Summed Poisson distribution as a function of At. 
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Since one is dealing with a discrete distribution function, the inequality 
is necessary and one usually obtains a value of x* that corresponds to a value 
slightly above P(0. 99) . This discrepancy is usually circumvented by the 
proper use of the words "at least" and "at most". For example, if At = 1, 
a probability table for exactly X events and the cumulative distribution is 
given in Table 5. From Table 5 one can see that for 98. 1 percent of the time 
there will be three or less events when At = i. 0, or for 1.9 percent of the time 
there will be four or more events. 

TABLE 5. POISSON DISTRIBUTION Next, the application of the 

FUNCTION (At = 1) above information to our problem of 

predicting the solar proton dose that 
would be seen on a space mission at 
various percentile levels will be ex- 
amined. For simplicity, assume that 
if an event occurs , it is depicted by the 
so-called "large" event model or model 
I of Table 3. Using the model and the 
sample statistical estimate of A^ = 0. 01 

event /week and for convenience assum- 
ing that the duration of exposure to 
large solar events is 100 weeks, then 
the value of At is 1.0 for the Poisson 
distribution ( Table 5) . Now assume 
that one wishes to find the dose levels 
that would be exceeded 27 percent, 2 percent, and 0. 1 percent of the time. 
Applying equation (25) with P taking on the values 0. 73, 0. 98, and 0. 999, 
it is seen from Table 5 that the values of x are 1, 3, and 5 events. There- 
fore, to find the dose at the respective levels, the doses behind the various 
shield thicknesses of model I are multiplied by the values of x above. The 
results are shown in Table 6 and plotted in Figure 3. Note that the 73. 6- 
percent level corresponds in our model to the mean or expected value of the 
Poisson distribution function. The percentile level at the expected value is 
found only when At is an integer. It becomes smaller, approaching 50 per- 
cent, as At increases to large values. For example, when At = 20, the per- 
centile at x* = 20 is 55. 9 percent [from equation (22)]. 

A similar computation can be carried out using event model II. The 
values of the Poisson mean become respectively A t = i. 0, A t = 2. 0, and 

Agt = 3. 0, where t = 100 weeks. Tables similar to Table 5 are constructed 

and the following results are found at the expected value, the 98-percent, and 


X 

P(x) 

2P(x) 

0 

0. 3679 

0. 3679 

l 

0. 3679 

0. 7358 

2 

0. 1839 

0. 9197 

3 

0.0613 

0. 9810 

4 

0. 0153 

0. 9963 

5 

0. 0031 

0. 9994 

6 

0. 0005 

0. 9999 
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SOLAR PROTON DOSE (RADS) (CENTER OF SPHERE) 



Figure 3. Solar proton dose curves as a function of shield thickness at 
various probability levels for 100-week missions. 
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TABLE 6. SOLAR PROTON DOSES (RADS) AT VARIOUS PROBABILITY 
LEVELS FOR 100-WEEK MISSIONS USING MODEL I 



x — 1 

x = 3 

x = 5 

Z ( g/ cm 2 ) 

7 3. 6-Percent 
Level 

98. 1 -Percent 
Level 

99. 9-Percent 
Level 

2 

1420 

4260 

7100 

5 

340 

1020 

1700 

10 

93 

279 

465 

15 

43 

129 

215 

20 

25 

75 

125 

30 

11 

33 

L_ 55 


the 99. 9-percent levels. For the large events, x = 1, 3, and 5; for the medium 
events, x = 2, 5, and 8; and for the small events, x = 3, 7, and 10. The per- 
centile values at the expected or mean value levels (x* = At) are 73. 6 per- 
cent, 67. 7 percent, and 64. 7 percent, respectively. To reduce the necessity 
of repetition, this set of levels shall be denoted as the approximately 70-per- 
cent level. The results are shown in Table 7. The comparison of results 
from the total columns of Table 7 with the comparable percentile level of 
Table 6 indicates very good agreement. The larger values of Table 7 at the 
70-percent level reflect the larger values of model II for 1 eventweek, shown 
in Table 3. However, at the 99. 9-percent level, model II seems to under- 
estimate the large model I event. This is because the five events of model I 
are larger than the five large events of model II, and the medium and small 
events of model II do not compensate in full for this difference. It is seen 
that for an exposure period of 100 weeks, model I is completely adequate for 
predicting the solar proton hazard at any reasonably high level of confidence 
(> 70 percent) . In general, the writer believes that if a level of confidence 
is chosen that leads to at least one large event, then model I should be suffi- 
cient for predicting the solar proton hazard. An example where this does 
not occur is when t is only 10 weeks; and if one calculates the 90-percent 
level, the value X = 0 satisfies the limits of equation (25) for the large event 
case (At = 0. 1) . For the medium (At = 0. 20) and small (At = 0. 30) event 
case, the value X = 1 satisfies the limits of equation (25). Thus, if model I 
is used there are no events q,t the 90-percent level; whereas, if model jH is 
used, one expects one medium event and one small event. Hence, for a shield 
of 5 g/cm 2 of aluminum, the point dose at the center of a spherical shell is 
0. 0 rads for model I and is 27. 4 rads for model II at the 90-percent level. 
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TABLE 7. SOLAR PROTON DOSES AT VARIOUS PROBABILITY LEVELS 
FOR 100-WEEK MISSIONS USING MODEL II 


a. Approximately 70-Percent Level 


Z(g/cm 2 ) 

x = 1 

x = 2 

x = 3 


D 1 

D 2 

D 3 

Total 

2 

1195 

190 

45 

1430 

5 

286 

48 

10. 2 

344. 2 

10 

78 

15 

2.7 

95. 7 

15 

36 

7.0 

1. 29 

44. 3 

20 

21 

4. 0 

0. 75 

25. 75 

30 

9. 3 

1. 6 

0. 33 

11. 23 


b. 98- 

Percent Level 


Z(g/cm 2 ) 

x = 3 

x = 5 

x = 7 


°1 

°2 

°3 

Total 

2 

3585 

475 

105 

4165 

5 

858 

120 

23. 8 

1002 

10 

234 

37. 5 

6. 3 

277. 8 

15 

108 

17. 5 

3. 01 

128. 5 

20 

63 

10 

1. 75 

74. 75 

30 

27. 9 

4 

0. 77 

32. 70 


c. 99. 9 -Percent Level 


Z (g/ cm 2 ) 

x = 5 

D i 

x = 8 

D 2 

X 

O II 

o 

Total 

2 

59 r 75 

760 

150 

6885 

5 

1430 

192 

34 

1656 

10 

390 

60 

9 

459 

15 

180 

2 a 

4. 3 

212 

20 

105 

16 

2. 5 

124 

30 

46. 5 

6. 4 

1. 1 

54 
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At the 99-percent level, one finds that model I predicts one large event, while 
model II predicts one large event, two medium events, and two small events. 
For the same shield as above , the model I dose is 340 rads and the model II 
dose is 341 rads, which is very good agreement at the 99-percent level. For 
the reasons given above, a useful rule would be: 

Empirical Rule : If at a given probability level there is at least one 

large event, model I is sufficient to describe the radiation hazard. 

If no large event is found at the given probability level, model II should 

be used. 

At this point it might be useful to depict several curves similar to 
Figure 3 for various mission times t, with A = 0. 01 event/week, at given prob- 
ability or percentile levels. To do this, arrays similar to Table 5 could be con- 
structed for various values of At, and the number of events, x, necessary to 
satisfy the inequality of equation (25) could be found for various values of P < 1. 

Rather than provide a multitude of similar graphs , Table 8 will allow 
the reader to find the number of events for a range of anticipated values of At 
at nine different percentile levels from 50 percent to 99. 9 percent. Table 8 
was constructed by choosing a probability level P, the number of events N, 
and then finding the value of m(At) that satisfied the following equation: 

N -m x 

Pr (x — N) = V — P- = P , (m = At) . (26) 

x =0 X - 

For example, if P = 0. 99 and N = 2, one has 


f(m) 


-m -m 

= e + me + 


mr 


-m 

e 


2 


0. 99 = 0 . 


(27) 


Equation (27) has the solution of m = 0. 44, which is found by Newton's method 
of successive iteration. Newton's method is simply, 


f (nik) 

m k+i = m k " f« (ny 


(28) 


The derivative of equation (26) can be shown to be simply, 


d 

dm 


l 


x=0 


-m x 
e m 


x : 



= f'(m) 


N -m 
-m e 

N! 


(29) 
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TABLE 8. THE VALUES OF Xt TO GIVE N OR LESS EVENTS 
AT THE EXACT PERCENTILES 


a. 

0. 001 

0. 005 

0. 010 

0. 025 

0.050 

0. 100 

0. 15 

0. 25 

0. 50 

^p a 

99.9 

99. 5 

99. 0 

97. 5 

95. 0 

90. 0 

85.0 

75. 0 

50. 0 


Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

Percent 

0 

0. 0010 

0. 0050 

0. 010 

0. 025 

0. 051 

0. 105 

0. 16 

0. 29 

0. 69 

1 

0. 045 

0. 10 

0. 15 

0. 24 

0. 36 

0. 53 

0. 68 

0. 96 

1.68 

2 

0. 19 

0. 34 

0. 44 

0. 62 

0. 82 

1. 10 

1. 33 

1. 73 

2. 67 

3 

0.43 

0. 67 

0. 82 

1.09 

1.37 

1.74 

2. 04 

2. 54 

3. 67 

4 

0.74 

1. 08 

1. 28 

1. 62 

1. 97 

2. 43 

2. 79 

3. 37 

4. 67 

5 

1. 11 

1. 54 

1. 79 

2. 20 

2. 61 

3. 15 

3. 56 

4. 22 

5. 67 

6 

1. 52 

2.04 

2. 33 

2. 81 

3. 29 

3. 89 

4. 35 

5. 08 

6. 67 

7 

1. 97 

2. 57 

2. 91 

3. 45 

3. 98 

4. 66 

5. 15 

5. 96 

7. 67 

8 

2. 45 

3. 13 

3.51 

4. 12 

4. 70 

5. 43 

5. 97 

6. 84 

8. 67 

9 

2. 96 

3.72 

4. 13 

4. 80 

5. 43 

6. 22 

6. 80 

7. 73 

9. 67 

10 

3. 49 

4. 32 

4. 77 

5. 49 

6. 17 

7. 02 

7. 64 

8. 62 

10. 67 

11 

4. 04 

4. 94 

5. 43 

6. 20 

6. 92 

7. 83 

8. 48 

9. 52 

11. 67 

12 

4. 61 

5. 58 

6. 10 

6. 92 

7. 69 

8. 65 

9. 34 

10. 42 

12. 67 

13 

5. 20 

6. 23 

6. 78 

7. 65 

8. 46 

9. 47 

10. 19 

11. 33 

13. 67 

14 

5. 79 

6. 89 

7. 48 

8. 40 

9. 25 

10. 30 

11. 06 

12. 24 

14. 67 

15 

6. 41 

7. 57 

8. 18 

9. 15 

10. 04 

11. 14 

11. 92 

13. 15 

15. 67 

16 

7. 03 

8. 25 

8. 89 

9. 90 

10. 83 

11. 98 

12.79 

14. 07 

16. 67 

17 

7. 66 

8. 94 

9. 62 

10. 67 

11. 63 

12. 82 

13. 67 

14. 99 

17. 67 

18 

8. 31 

9. 64 

10. 35 

11.44 

12. 44 

13. 67 

14. 55 

15. 91 

18. 67 

19 

8. 96 

10. 35 

11. 08 

12. 22 

13. 25 

14. 53 

15. 43 

16. 83 

19. 67 

20 

9. 62 

11. 07 

11. 83 

13. 00 

14. 07 

15. 38 

16. 31 

17.75 

20. 67 

21 

10. 29 

11.79 

12. 57 

13.79 

14. 89 

16. 24 

17. 20 

18. 68 

21. 67 

22 

10. 96 

12. 52 

13. 33 

14. 58 

15. 72 

17. 11 

18. 09 

19. 61 

22. 67 

23 

11. 65 

13. 26 

14. 09 

15. 38 

16. 55 

17. 97 

18. 98 

20. 54 

23. 67 

24 

12. 34 

14. 00 

14. 85 

16. 18 

17. 38 

18. 84 

19. 88 

21. 47 

24. 67 

25 

13. 03 

14.74 

15. 62 

16.98 

18. 22 

19. 72 

20. 77 

22. 40 

25. 67 

26 

13. 73 

15. 49 

16. 40 

17. 79 

19. 06 

20. 59 

21. 67 

23. 34 

26. 67 

27 

14. 44 

16. 25 

17. 17 

18. 61 

19. 90 

21. 47 

22. 57 

24. 27 

27. 67 

28 

15. 15 

17. 00 

17. 96 

19. 42 

20. 75 

22. 35 

23. 48 

25. 21 

28. 67 

29 

15. 87 

17.77 

18.74 

20. 24 

21. 59 

23. 23 

24. 38 

26. 15 

29. 67 

30 

16. 59 

18. 53 

19. 53 

21. 06 

22. 44 

24. 11 

25. 28 

27. 09 

30. 67 

31 

17. 32 

19. 30 

20. 32 

21. 89 

23. 30 

25. 00 

26. 19 

28. 03 

31.67 

32 

18.05 

20. 08 

21. 12 

22.72 

24. 15 

25. 89 

27. 10 

28. 97 

32. 67 

33 

18.78 

20. 86 

21. 92 

23. 55 

25. 01 

26. 77 

28. 01 

29. 91 

33. 67 

34 

19. 52 

21. 64 

22.72 

24. 38 

25. 87 

27. 66 

28. 92 

30. 85 

34. 67 

35 

20. 26 

22. 42 

23.53 

25. 21 

26.73 

28. 56 

29. 83 

31. 79 

35. 67 

36 

21. 00 

23. 21 

24. 33 

26. 05 

27. 59 

29. 45 

30.75 

32. 74 

36. 67 

37 

21.75 

24. 00 

25. 14 

26. 89 

28. 46 

30. 34 

31. 66 

33. 68 

37. 67 

38 

22. 51 

24.79 

25. 96 

27.73 

29. 33 

31. 24 

32. 58 

34. 63 

38. 67 

39 

23. 26 

25. 59 

26. 77 

28.58 

30. 20 

32. 14 

33. 50 

35. 57 

39. 67 

40 

24. 02 

26. 38 

27. 59 

29. 42 

31.07 

33. 04 

34. 42 

36. 52 

40. 67 


a. P = 1 - a 

b. If Xt lies between two entries, use N corresponding to the largest entry. Thus, if 
Xt = 0. 5, use N = 4 at the 99. 9-percent level. If Xt is less than or equal to the top 
entry, use N = 0. Thus, if Xt = 0. 02, N ■= 0 at the 97. 5-percent level, but N = 1 at 
the 99. 0-percent level. 
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From equations (26) and (29), the iteration becomes 


“k +i =m k + 



( 30 ) 


The value of m = At found in this manner is the expected value or mean of the 
Poisson distribution that has a percentile level of exactly P x 100 percent for 
N or less events. For example, if P = 0. 99, and N = 2, At = 0. 44 satisfies 
equation (26) . It should be noted that if At = 0. 70, N = 3 must be chosen at the 
99-percent level. This follows from equation (25) , and P is no longer exactly 
99 percent but becomes 99. 42 percent. Hence, one says that for At = 0. 70, the 
probability is at least 0. 99 that the number of events is three or less. Another 
way of stating the condition is that the probability is no greater than ot = 0. 01 
that there will be more than three events when At = 0. 70; i. e. , 1 - P = 0. 0058. 


The top of Table 8 is headed with a row of values labeled ot , which de- 
notes the probability of more than N events or simply ot = 1 - P. From 
equation (23) it can be seen that Pr(X>N) = 1 - P(X^N) = ot. The probability 
statements using ot will be made as follows: the probability is no greater 
than ot that more than N events will be observed in t weeks. The use of 
the ot -probabilities will be derived in the next section. The range of values 
in Table 8 should provide for the refinements necessary in the following work. 

If one uses the discrete Poisson distribution properly, it is not possible 
to get better values for N events at the preassigned probability levels than 
those obtained with the prescription used for Table 8. It should be remembered 
that a discrete distribution cannot have fractional events. The user should 
not be tempted to interpolate between probability values of the cumulative dis- 
tribution function and obtain a fractional value of X that is used in the same 
manner as the writer is using discrete events in order to meet some special 
requirement. There are certain practical exceptions to this rule , but they 
occur only in that domain where the Poisson is well approximated by the 
Gaussian or normal distribution. This does not occur until At > 20. 


It might be of value to examine the case where At is sufficiently large 
that the Gaussian distribution becomes a reasonable approximation to the 
Poisson distribution. Where this becomes a safe assumption depends on the 
user’s threshold of acceptable accuracy. For example, even at a value of 
m = 40, the probability of having a value of X ^ 40 is still 54 percent, where- 
as for a Gaussian distribution, the area to the left of the mean is exactly 
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50 percent. The Poisson distribution always has a slight skew to the right even 
for very large values of m = At. However, with symmetric distributions such 
as the Gaussian, it is common to use measures of central tendency about the 
mean for making probability statements. If the Poisson is to be approximated 
by a normal distribution, the value of the mean is m = At and the standard 
deviation is sfm for the "equivalent" Gaussian distribution. The transforma- 
tion to the normal distribution in the variable t is given by 

i _tV 2 

0(t)dt = e / dt , (31) 

sf2 n 


where 

x - m 

l • 

si m 

Here, the mean value of the normal distribution is at t = 0, and the standard 
deviation cr = 1, x is the number of events for the Poisson, and m is the 

l 

mean of the Poisson (m = At). The values of the integral, 
t 

/ 0(t)dt (32) 

-t 

are available in most handbooks , where t = 1 is one standard deviation and 
multiples of t are thought of as multiples of the standard deviation. 

The reader will find that if probabilities are measured about the mean 
rather than from the tails, the Poisson distribution is fairly well approximated 
by the Gaussian even down to relatively low values of At (e. g. , 25). (See 
Figure 6 of Section IV for a heuristic comparison. ) 

It is of interest to recall that in dealing with a discrete distribution 
function, the difference between two cumulative distribution or sums [for 
example, F(X") and F(X’)J denotes the probability: 

P(X' < X ^ X”) = F(X") - F(X') . (33) 

As an example, the probability f(x) of obtaining x dots when two dice are 
rolled is given in Table 9. 
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TABLE 9. PROBABILITIES FOR DICE PAIRS 


X 

2 3 4 5 

6 

7 

8 

9 

10 

11 

12 

f(x) 

1/36 2/36 3/36 4/36 

5/36 

6/36 

5/36 

4/36 

3/36 

2/36 

1/36 

F(x) 

1/36 3/36 6/36 10/36 

15/36 

21/36 

26/36 30/36 33/36 

35/36 

36/36 


Therefore, from equation (33) , F(7) - F(4) = P(4 < X s 7 ) = 15/36, or the 
probability of getting a 5, 6 , or 7 on a pair of dice is 15/36. 


SECTION IV. CONFIDENCE INTERVALS ARE ESTABLISHED 


The results of the previous section depended strongly on the choice 
of the parameter A or the probability of having a large event week. In this 
section the confidence that can be placed on the value of A as calculated from 
observed sample data will be established. One of the virtues of the Poisson 
distribution is that the value of the mean ( m = At) completely determines the 
distribution function (i. e. , cr 2 = m) ; whereas the mean and variance (a 2 ) are 
needed for most distribution functions, and they are not simply related. Thus, 
in working with small sampling statistics from a normal distribution, one 
needs to establish confidence intervals for both the mean and the variance. 

To more clearly explain what is being done in this section, an example 
will be given. The Poisson law arises very often in certain biological prob- 
lems, such as with organisms distributed at random over the bottom of a lake. 
The number of such organisms found in a series of trial dredgings from sepa- 
rate small areas of the same size will follow this law. Statisticians calculate 
boundaries of possible outcomes from a given small sample, and these values 
are called confidence limits at a certain probability level for the assumed 
distribution function. If a biologist counted 2i organisms from one of his 
dredgings , he could assert that he is 95-percent confident that the mean or 
expected value lies between 13 and 32 organisms per unit area assuming a 
Poisson distribution. Thus, with only one sample and the assumption of a 
Poisson process, it is possible to set upper and lower bounds on possible out- 
comes at a given probability level. From the above example, it can be asserted 
that if many dredgings are made , it is expected that only 5 percent will con- 
tain a number of organisms outside the predicted range. 

To be more specific for the class of problems herein, it will be 
assumed that the number of events observed over the sample space of T 
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weeks is given by X events. First, it should be made clear that a small 

A 

A 0 

sample estimate of A = exists. In the previous section it was assumed 

that m was the true mean for the Poisson distribution and that probability 
statements could be made;, 

Pr(X > , (34) 


where X represents the total number of events that may occur in a time of 

T weeks, a is some small risk such as 0. Of, or 0. 001, and N is the 

a 

smallest value of X that satisfies the inequality at the small risk. An upper 
bound on m (At) is needed now to find a value so that the statement, 

[pr(x>N^^o;J , (35) 

would be correct with a high probability; for example, 0. 99 or, in general, P. 


Thus, one finds a one-sided upper confidence interval pn the true value 
of m = At. Suppose the probability is P that one would see more than X^ 

events for the T weeks of observation. Then the true value of m would have 
to be a value such that 


Pr( x £ Xq 


) 


0 

2 

x=0 


-m x 
e m 


?£ ! 


1 - P 


(36) 


where i - P is a small number. 

Let m' be the solution to this equation. The will be that small- 
est integer for which 

Pr ^x > | m = m' ^ =£ a . (37) 

That is, is the smallest integer satisfying 


N 

a 

Z 

x=0 


-m' 


(m')' 


Xj 


a 


(38) 
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For example, if X Q = 3, P= 0. 99, and ol ~ 0. 001, then m' is the solution to 
3 -m' . x 

V _® L22J — = i _ 0. 99 = 0. 01 , (39) 

x=0 X * 


or 


-m 

e 



i + m 1 + 


( m ») 2 

2 


+ 


( m ») 3 
6 . 


0. 01 


(40) 


Using a table of the cumulative Poisson distribution or Newton's method of 
finding the root of a transcendental equation, one finds that m' = 10. 05 satis- 
fies the equation. This means that one is 99-percent certain (confident) that 
the true mean for large events in the 300 weeks is 10 or less. 

Next, one finds the smallest integer N^ that satisfies 
» 

N 

2: 1 - a = 0. 999 . (41) 


a 

Z 

x=0 


,-io 


10 


x , 


Again, by using tables of Poisson cumulative probabilities or Table 8, one 

finds that N f = N* = 21 is the solution. Some other values for N* when 

a 0. 001 ol 

P = 0. 99 (m T = 10) are readily found from Table 8 and are shown in Table 10. 


TABLE 10. SUMMARY OF VALUES N* 

OL 


P = 0. 99 

a 

0. 001 

0. 01 

0. 025 

0. 05 

0. 10 

m* = 40. 05 

N* 

a 

21 

18 

17 

15 

14 


As a summary of the above, one is 100P percent confident that the 
following statements are true 4 : 

The true value of m is no greater than m', and the probability 

is no greater than a that more than N^ event weeks will be observed 

in T weeks. 


4. Most of these concepts grew out of a private communication with Dr. James 
E. Norman, Jr. , Department of Statistics, University of Georgia. 
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Returning to the basic problem, one wishes to establish the confidence 
interval on m , hence to find the probability 

Pr (m" ^ m ^ m') ^ (1 - 2/3) 100% , (42) 

where 1 - P = /3. If P = 0. 99 and /3 = 0. 01, then 

Pr (m" ^ m ^ m') ^ 98% 


and one is at least 98 percent certain that m lies between m' and m" , where 
m has only a 1 -percent chance of being greater than m'. The general set of 
equations can now be written that will determine the value of m' and m" : 


0 

2 

x=0 



i - P = /3 


(43) 


and 


2 

x=X„ 


-m' 


(m") 


X 


= p 


or 


V 1 


2 

x=0 


-m' 


(m'T 


x.' 


1-/3 


(44) 


where is the observed number of events. The form of the above equation 

is an outgrowth of using a discrete distribution function and follows the com- 
mon practice of textbooks on statistics. Solving equations (43) and (44) for 
m* and m" by using Newton's method [see equation (30)] , the upper and 
lower bounds on m are found for X^ observed events. Table ii summarizes 

these results for eight different values of /3 or P. Figures 4 and 5 depict 
typical results taken from Table 11. 


Figure 6 presents comparable results when the normal law is assumed. 
For this calculation, the change of variable to the normal variate t is made 
using the relationship, 


t 


X„ - x 


J 


X, 


(45) 
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TABLE il. UPPER AND LOWER BOUNDS FOR THE MEAN 
AT PROBABILITY OF P = (i - 0) 100% 


Confidence 

Interval 

99. 8 
Percent 

99. 0 
Percent 

98.0 

Percent 

95 

Per 

0. ( 

. 0 

cent 

P 

0. 001 

0. 005 

0. 010 

)25 

Events, 

Lower 

Bound a 

Upper 

Bound 

Lower 

Bound 

Upper 

Bound 

Lower 

Bound 

Upper 

Bound 

Lower 

Bound 

Upper 

Bound 

0 

0 

6. 91 

0 

5. 30 

0 

4. 61 

0 

3.69 

i 

0. 001 

9. 23 

0.005 

7. 43 

0. OK 

6. 64 

0. 025 

5. 57 

2 

0. 045 

11. 23 

0. 10 

9. 27 

0. 15 

8. 41 

0. 24 

7. 22 

3 

0. 19 

13. 06 

0. 34 

10. 98 

0. 44 

10. 05 

0.62 

8. 77 

4 

0. 43 

14. 79 

0. 67 

12. 59 

0. 82 

11. 60 

1. 09 

10. 24 

5 

0. 74 

16. 45 

1.08 

14. 15 

1. 28 

13. 11 

1. 62 

11. 67 

6 

1. 11 

18. 06 

1. 54 

15. 66 

1. 79 

14. 57 

2. 20 

13. 06 

7 

1. 52 

19. 63 

2. 04 

17. 13 

2. 33 

16. 00 

2. 81 

14.42 

8 

1. 97 

21. 16 

2. 57 

18. 58 

2. 57 

17. 40 

3. 45 

15.76 

9 

2. 45 

22. 66 

3. 13 

20. 00 

3. 51 

18. 78 

4. 12 

17. 08 

10 

2. 96 

24. 13 

3.72 

21. 40 

4. 13 

20. 14 

4. 80 

18. 39 

11 

3. 49 

25. 59 

4. 32 

22. 78 

4. 77 

21. 49 

5. 49 

20. 96 

12 

4. 04 

27. 03 

4. 94 

24. 15 

5. 43 

22. 82 

6. 20 

20. 96 

13 

4. 61 

28. 45 

5. 58 

25. 50 

6. 10 

24. 14 

6. 92 

22. 23 

14 

5. 20 

29. 85 

6. 23 

26. 84 

6. 78 

25. 45 

7. 65 

23. 49 

15 

5. 79 

31. 24 

6. 89 

28. 16 

7. 48 

26. 74 

8. 40 

24.74 

16 

6. 41 

32. 62 

7. 57 

29. 48 

8. 18 

28. 03 

9. 15 

25. 98 

17 

7. 03 

33. 99 

8. 25 

30. 79 

8. 89 

29. 31 

9. 90 

27. 22 

18 

7. 66 

35. 35 

8. 94 

32. 09 

9. 62 

30. 58 

10. 67 

28. 45 

19 

8. 31 

36. 70 

9. 64 

33. 38 

10. 35 

31. 85 

11.44 

29. 67 

20 

8. 96 

38. 04 

10. 35 

34. 67 

11. 08 

33. 10 

12. 22 

30. 89 

21 

9. 62 

39. 37 

11. 07 

35. 95 

11. 83 

34. 35 

13. 00 

32. 10 

22 

10. 29 

40. 70 

11.79 

37. 22 

12. 57 

35. 60 

13. 79 

33. 31 

23 

10. 96 

42. 02 

12. 52 

38. 48 

13. 33 

36. 84 

14. 58 

34. 51 

24 

11. 65 

43. 33 

13. 26 

39. 74 

14. 09 

38. 08 

15. 38 

35. 71 

25 

12. 34 

44. 64 

14.00 

41. 00 

14. 85 

39. 31 

16. 18 

36. 90 


a. The lower bound is calculated so that j 8 is less than the indicated 
value. 
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TABLE 11. (Concluded) 


Confidence 

90. 

0 

80. 

0 

70 

.0 

50 

.0 

Interval 

Percent 

Percent 

Percent 

Percent 

£ 

0. 050 

0. 100 

0. 

15 

0. 

25 

Events, X Q 

Lower 

Upper 

Lower 

Upper 

Lower Upper 

Lower Upper 

Bound 

Bound 

Bound 

Bound 

Bound 

Bound 

Bound 

Bound 

0 

0 

3. 00 

0 

2. 30 

0 

1. 90 

0 

1. 39 

1 

0. 051 

4. 74 

0. 105 

3. 89 

0. 16 

3. 37 

0. 29 

2. 69 

2 

0. 36 

6. 30 

0. 53 

5. 32 

0. 68 

4. 72 

0. 96 

3. 92 

3 

0. 82 

7.75 

1. 10 

6. 68 

1. 33 

6. 01 

1. 73 

5. 11 

4 

1. 37 

9. 15 

1.74 

7. 99 

2. 04 

7. 27 

2. 54 

6. 27 

5 

1. 97 

10. 51 

2. 43 

9. 27 

2. 79 

8. 49 

3. 37 

7. 42 

6 

2. 61 

11. 84 

3. 15 

10. 53 

3. 56 

9. 70 

4. 22 

8. 56 

7 

3. 29 

13. 15 

3. 89 

11. 77 

4. 35 

10. 90 

5. 08 

9. 68 

8 

3. 98 

14. 43 

4. 66 

12. 99 

5. 15 

12. 08 

5. 96 

10. 80 

9 

4. 70 

15. 71 

5. 43 

14. 21 

5. 97 

13. 25 

6. 84 

11. 91 

10 

5. 43 

16. 96 

6. 22 

15. 41 

6. 80 

14. 41 

7. 72 

13. 02 

11 

6. 17 

18. 21 

7. 02 

16. 60 

7. 64 

15. 57 

8. 62 

14. 12 

12 

6. 92 

19. 44 

7. 83 

17. 78 

8. 48 

16. 71 

9. 52 

15. 22 

13 

7. 69 

20. 67 

8. 65 

18. 96 

9. 34 

17. 86 

10. 42 

16. 31 

14 

8. 46 

21. 89 

9. 47 

20. 13 

10. 19 

19. 00 

11. 33 

17. 40 

15 

9. 25 

23. 10 

10. 30 

21. 29 

11. 06 

20. 13 

12. 24 

18. 49 

16 

10. 04 

24. 30 

11. 14 

22. 45 

11. 92 

21. 26 

13. 15 

19. 57 

17 

10. 83 

25. 50 

11. 98 

23. 61 

12. 79 

22. 38 

14. 07 

20. 65 

18 

11. 63 

26. 69 

12. 82 

24. 76 

13. 67 

23. 50 

14. 99 

21. 73 

19 

12. 44 

27. 88 

13. 67 

25. 90 

14. 55 

24. 62 

15. 91 

22. 81 

20 

13. 25 

29. 06 

14. 53 

27. 05 

15. 43 

25.74 

16. 83 

23. 88 

21 

14. 07 

30. 24 

15. 38 

28. 18 

16. 31 

26. 85 

17.75 

24. 96 

22 

14. 89 

31. 41 

16. 24 

29. 32 

17. 20 

27. 96 

18. 68 

26.03 

23 

15. 72 

32. 59 

17. 11 

30. 45 

18. 09 

29. 07 

19. 61 

27. 10 

24 

16. 55 

33. 75 

17. 97 

31. 58 

18. 98 

30. 17 

20. 54 

28. 17 

25 

17. 38 

34. 92 

18. 84 

32.71 

19. 88 

31. 28 

21. 47 

29. 23 
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Figure 4. 99. 8-percent confidence limits for Poisson law. 
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Figure 5. 90. 0-percent confidence limits for Poisson law. 
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Figure 6. Comparison of 98. 0-percent confidence limit 
for the Poisson and normal laws. 

where is the observed number of events. (See equation ( 31) for an ex- 
planation. ) 

Even though confidence bounds on the estimate of m = At have been ob - 
tained, the major interest is with the upper bounds that an investigator should 
use when only events have been observed over a T week time interval. 

Following this concept and the simple model of a large dose event week, it is 
recalled that only three such event weeks occurred over the 300 active weeks 
of the 19th solar cycle. Now, if one wishes to find the upper bounds on this 
observation, he may use Table 11, in conjunction with Table 8, and arrive 
at the lOOP-percent confidence level that the probability that events will 

occur is no greater than some small value a. Thus, Table 12 is constructed. 
This table is a detailed extension of the results given in Table 10 with different 
values of P and m* for different values of a and N . 

The results shown in Table 12 are fairly clear; however, it may be use- 
ful to explain how Table 12 was constructed. First, the number of observed 
events X Q = 3 is used as the starting point. Using Table 11 suppose the 

95-percent upper bound on m is desired, then the value of m' = 7. 75 is 
found opposite X^ = 3 in the column /3 = 0. 05. Next, one refers to Table 8 
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and under the column headed by a = 0. 050, the value of m' = At = 7. 75 is found 
between 7. 69 and 8. 46 that corresponds to the N-numbers (first column) of 

12 and 13. Since the value of N is discrete, one must choose N = 13 to 

a 

guarantee that the probability is less than 0. 05 that more than 13 events will 
occur. It should be clear to the reader that the exact value of a lies between 
0. 05 and 0. 025, as can be seen from Table 8. However, because of the rather 
extensive computations necessary, it is sufficient to say that the probability 
is no greater than a (0. 05) that more than (13) events will be observed in 
300 solar active weeks. 

Since the 20th solar cycle has produced no large events, it is of interest 
to ask. If zero large events occur for a period of 300 active weeks, what are 
the possible upper bounds of m' and the number of events that could be 

exceeded at some small probability? (This may be a nonsensical question and 
will be discussed in the next section. ) Table 13 is constructed with this in 
mind. The methods used are the same as for Table 12. It is of interest to 
note that even though X Q = 0 events are observed, the 95-percent upper bound 

for the Poisson mean is m* = 3, which was the actual observed number of large 
event weeks during the 19th cycle. Thus, it seems that for a reasonable level 
of confidence (95 percent) , the analyst would be justified if he used m' = 3 
(A = 0. 01 event/week) for large events even though the solar activity for a 
given cycle was considerably different from that of the 19th cycle. Perhaps, 
the only conclusion that can be drawn is that the observed 1 9th cycle dose 
events could be used for any near average solar cycle, and if a cycle is pre- 
dicted to be similar to the 19th cycle, the results of Table 12 (X^ = 3) should 
be seriously considered as a possible model. 

To expand the events observed to cover both medium and small events 
as described in model II (Table 3) , Tables 14 and 15 are presented for X Q = 6 
and 9 observed events. Finally, Table 16 is presented as a summary of 
values to use for A ( events/week) at various confidence levels. For large 
events , it seems that the X^ = 0 column is probably reasonable to use if the 

solar cycle is not very active. If a cycle similar to the 19th is forecast, then 
the column under X Q = 3 is preferred. 

The value of A and the dose event model to use during the remaining 
270 less-active weeks of a solar cycle are not obtainable from the present 
analysis. However, until more data are available, the best one can do is use 
some criteria that are essentially arbitrary. For example, if one believes 
that the chance of a large dose event week is definitely dependent on some 
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TABLE 12. THE NUMBER OF EVENTS N^ THAT WILL BE EXCEEDED 


AT A PROBABILITY NO GREATER THAN a FOR THE iOQ^PERCENT 
UPPER BOUND VALUE OF THE POISSON MEAN, m', 

WHEN X = 3 OBSERVED EVENTS 


f 

m 

100P 

Percent 





a 





0.001 

0.005 

0.010 

0.025 

0.050 

0.100 

0. 150 

0.250 

0.500 

13.06 

99.9 

25 

23 

22 

20 

19 

18 

17 

15 

13 

10.98 

99.5 

22 

20 

19 

18 

17 

15 

14 

13 

11 

10.05 

99.0 

21 

19 

18 

17 

15 

14 

13 

12 

10 

8.77 

97.5 

19 

17 

16 

15 

14 

13 

12 

11 

9 

7.75 

95.0 

18 

16 

15 

14 

I3 a 

11 

11 

9 

8 

6.68 

90.0 

16 

14 

13 

12 

11 

10 

9 

8 

7 

6. 01 

85.0 

15 

13 

12 

11 

10 

9 

8 

7 

6 

5. 11 

75.0 

13 

12 

11 

10 

9 

8 

7 

6 

5 

3.67 

50.0 

11 

9 

9 

8 

7 

6 

6 

5 

3 


a. Example: One is 95 percent confident that the true mean does not exceed 
7. 75 (X(f=3) , and using this value of the mean, one is certain that the 
probability of seeing more than 13 events during a 300 week active period 
is no greater than 0. 050. 


TABLE 13. THE NUMBER OF EVENTS n' THAT WILL BE EXCEEDED AT A 

a 


PROBABILITY NO GREATER THAN a FOR THE 100P-PERCENT UPPER 
BOUND VALUE OF THE POISSON MEAN, m', WHEN X Q = 0 
OBSERVED EVENTS 


! 

m 

100P 

Percent 

a 

0.001 

0. 005 

0. 010 

0.025 

0.050 

0. 100 

0. 150 

0.250 

6. 91 

99. 9 


15 

14 

12 

11 


10 

9 

5.30 

99.5 

IKS 

12 

11 


9 

8 

8 

7 

4.61 

99.0 

B 

11 


9 

8 

7 

7 

6 

3.69 

97.5 

B 

9 

9 

8 

7 

6 

6 

5 

3.00 

95.0 

H 

8 

8 

7 

6 

5 

5 

4 

2.30 

90.0 

8 

7 

6 

6 

5 

4 

4 

3 

1.90 

85.0 

7 

6 

6 


4 

4 

3 

3 

1.39 

75.0 

6 

5 



4 

3 

3 

2 

0.69 

50.0 

4 

4 



2 

2 

2 

1 


0.500 

7 

5 


3 

2 

2 

1 

0 
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TABLE 14. THE NUMBER OF EVENTS N' THAT WILL BE EXCEEDED AT A 


PROBABILITY NO GREATER THAN a FOR THE 100P-PERCENT UPPER 
BOUND VALUE OF THE POISSON MEAN, m', WHEN X =6 
OBSERVED EVENTS 0 


t 

m 

100P 

Percent 





a 





0.001 

0.005 

0.010 

0.025 

0.050 

0.100 

0. 150 

0.250 

0.500 

18.06 

99.9 

33 

30 

29 

27 

25 

24 

22 

21 

18 

15.66 

99.5 

29 

27 

26 

24 

22 

21 

20 

18 

15 

14.57 

99.0 

28 

25 

24 

22 

21 

20 

19 

17 

14 

13.06 

97.5 

26 

23 

22 

21 

19 

18 

17 

15 

13 

11.84 

95.0 

24 

22 

21 

19 

18 

16 

15 

14 

12 

10.53 

90.0 

22 

20 

19 

17 

16 

15 

14 

13 

10 

9.70 

85.0 

21 

19 

18 

16 

15 

14 

13 

12 

10 

8. 56 

75.0 

19 

17 

16 

15 

14 

12 

12 

10 

8 

6.67 

. 

50.0 

16 

14 

13 

12 

11 

10 

9 

8 

6 


TABLE 15. THE NUMBER OF EVENTS N' THAT WILL BE EXCEEDED AT A 

a 

PROBABILITY NO GREATER THAN a FOR THE 100P-PERCENT UPPER 
BOUND VALUE OF THE POISSON MEAN, m’, WHEN X =9 
OBSERVED EVENTS 0 



100P 

Percent 


1 

i 

L. 

| a 

l 

i 



T 

m 

0.001 

0.005 

0.010 

0.025 

0.050 

0. 100 

0. 150 

0.250 

0.500 

22.66 

99. 9 

39 

36 

34 

32 

31 

29 

28 

26 

22 

20.00 

99.5 

35 

32 

31 

29 

28 

26 

25 

23 

20 

18.78 

99.0 

33 

31 

30 

28 

26 

24 

23 

22 

19 

17.08 

97.5 

31 

29 

27 

26 

24 

22 

21 

20 

17 

15.71 

95.0 

29 

27 

26 

24 

22 

21 

20 

18 

16 

14.21 

90.0 

27 

25 

24 

22 

21 

19 

18 

17 

14 

13.25 

85.0 

26 

23 

22 

21 

19 

18 

17 

16 

13 

11.91 

75.0 

24 

22 

21 

19 

18 

16 

15 

14 

12 

9.67 

50.0 

21 

19 

18 

16 

15 

14 

13 

12 

9 


M439 


39 



TABLE 16. VALUES FOR A ? = m'/300 (EVENTS/WEEK) FOR X 
OBSERVED EVENTS AT THE 100P-PERCENT 
UPPER BOUND CONFIDENCE LEVEL 


Vo 

100P\ 

Percent^ 

0 

3 

6 

9 

99. 9 

0. 0230 

0. 0435 

0. 0602 

0. 0765 

99. 5 

0. 0177 

0. 0366 

0. 0522 

0. 0667 

99. 0 

0. 0154 

0. 0335 

0. 0486 

0. 0626 

97. 5 

0. 0123 

0. 0292 

0. 0435 

0. 0569 

95. 0 

0. 0100 

0. 0258 

0. 0395 

0. 0524 

90. 0 

0. 0070 

0. 0223 

0. 0351 

0. 0474 

85. 0 

0. 0063 

0. 0200 

0. 0323 

0. 0442 

75. 0 

0. 0046 

0. 0170 

0. 0285 

0. 0392 

50. 0 

0. 0023 

0. 0122 

0. 0222 

0. 0322 

OBS a 

0. 0000 

0. 0100 

0. 0200 

0. 0300 


a. OBS denotes the value of X^/300. 

minimal level of solar activity and this level is not approached during the quiet 
periods of the cycle, it would be unreasonable to use the large event dose week 
even at a very low probability level , such as the values shown under the X^ = 0 
column of Table 16. Even if A is only 0. 002 event/week for the Poisson 
distribution, one sees that the chance of getting two or more large events is 
approximately 6 percent for T = 200 weeks. However, one cannot preclude 
the possibility of some smaller dose event occurring for the quiet period of a 
solar cycle. This writer suggests using A = 0. 02 event/week for the quiet 
period, but recommends that the dose model for the small dose event (D^) of 
model n (Table 3) be used. This is equivalent to expecting about one such 
event per year, which is reasonably close to the actual observed number of 
events from October 1961 through July 1966. 

To illustrate applications of the foregoing work, eight different trip 
lengths will be considered ranging from 13 weeks to 260 weeks (5 years) , 
during the most active 6 years of a solar cycle that is similar to the 19th 
cycle (X Q = 3). To simplify the possible combinations, four different prob- 
ability levels for the upper bound mean (At) are used at four probability levels 
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for the Poisson distribution. A mean (0. Oit) corresponding to the actual ob- 
served (OBS) large event weeks of the 19th cycle is also given. This corre- 
sponds to the 95-percent level of the mean for = 0 events, and this value is 

recommended for the average type cycle. The results of these combinations 
are shown in Table 17. The entries in Table 17 give the number of large events 
that will be exceeded at a probability equal to or less than a, corresponding 
to the Poisson mean (At) that one is 100P percent certain will not be exceeded 
(see the footnote for Table 12). Table 17 is constructed in a manner similar to 
Table 12; however, one begins with Table 16 for values of A at various levels 
of percent of confidence (100P). 

The foregoing method is very convenient, and one can readily generate 
many such tables. However, it would be instructive to demonstrate to the 
reader that any more precise approach would be excessively tedious and the 
probabilities would be different for each combination of m(At) and the number 
of events N. By using a Poisson distribution table, one can arrive at the 
exact probabilities of exceeding N events for a given value of m = At. These 
exact values were calculated for a 78-week mission using the 95-percent value 
of m(At = 2. 00) . These exact results are shown in parentheses in Figure 7. 
The first quantities given, which identify the percent probabilities of Figure 7, 
are taken from the a values of Table 17. It should be noted that the methods 
used in forming Table 17 give a conservative estimate of the true probabilities 
shown in parentheses. 

Figures 8, 9, and 10 illustrate values of a at 0. 001, 0. 01, and 0. 1 
for different upper bound values of the mean showing a range from the ob- 
served (approximately 50 percent) to the 99. 9-percent confidence level of 
the mean (At) for a 78-week mission during the active weeks of a very active 
cycle. Comparisons are also shown in Figures 8, 9, and 10 to the work of 
other authors [2, 9, 10] who have made similar computations but have used 
different models for predicting the solar proton dose. An interesting aspect 
of the above comparison is that at the 0. 1 -percent (100 x a) level, the pres- 
ent work is considerably lower even for the 99. 9-percent upper bound value 
of At. However, at the 10-percent level, the reverse situation seems to be 
the case. The major difference in the methods of most of the other authors 
and the present work is that the size of our large event is fixed but the number 
of events may be quite large , whereas in the other methods the size of a single 
event may be extremely large (several times larger than the large event used 
in this work) . These differences will be discussed in the last section. 

Figure 11 illustrates the conjectured proton event dose at the center 
of a spherical shell of aluminum during the quiet periods of a solar cycle. 

In this calculation, A = 0. 02 event/week and the product At is found in the 
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TABLE 17. NUMBER OF LARGE EVENT WEEKS EXPECTED FOR 
VARIOUS MISSION LENGTHS AT GIVEN JOINT PROBABILITIES 

P AND a FOR X = 3 


a. 

Xt 

N 

Px 10C 

>7 

99.9 

14 

99.0 

U 

95.0 

19 

90. 0 

L3 

OBS 

b. 


a. t = 13 Weeks 


0.001 0. 01 0.05 0. 10 


e. t = 104 Weeks 


m = Xt 

a 

Px ioX 

0. 001 

0.01 

0. 05 0. 10 

4. 52 

99.9 

12 

10 

8 7 

3. 48 

99. 0 

10 

8 

7 6 

2. 68 

95.0 

9 

7 

6 5 

2. 32 

90.0 

9 

6 

5 4 

1* 04 

OBS 

5 

4 

3 2 


b. t = 26 Weeks 


156 Weeks 


Xt P x 100N 0. 001 0. 01 0. 05 0. 10 


3 2 

3 2 

2 1 


c. t = 52 Weeks 


Xt Pxl00\J 0.001 0.01 0.05 0.10 


6 99. 9 

4 99. 0 

4 95. 0 

6 90. 0 

2 OBS 


3 3 

3 3 

2 1 


m = Xt 

N \0! 
P x 100\ 

0.001 

0. 01 

0. 05 

0. 10 

6. 79 

99. 9 

16 

16 

11 

10 

5. 23 

99. 0 

14 

11 

9 

8 

4. 02 

95. 0 

11 

9 

8 

7 

3. 48 

90. 0 

10 

8 

7 

6 

1. 56 

OBS 

7 

5 

4 

3 

g. t 

= 208 Weeks 


d. t = 78 Weeks 


m = Xt 

\a 
Px loX 

0.001 

0. 01 

0. 05 , 

0. 10 

9.05 j 

99. 9 

20 

17 

14 

13 

6. 97 

99. 0 

16 

14 

12 

10 

5. 37 

95. 0 

14 

11 

9 

8 

4. 64 

90.0 

13 

10 

8 

7 

2.08 

OBS 

8 

6 

5 

4 


h. t 

= 260 Weeks 




\t I P x 100N 0. 001 0. 01 0. 05 0. 10 I 


5 5 

5 4 

4 3 

2 2 



At Px lOONJ 0. 001 0. 01 0. 05 0. 10 
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Figure 7. Solar proton dose curves as a function of shield thickness 
at various probability levels for a 78-week mission 
using a mean at the 95-percent level. 
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Figure 8. Solar proton dose curves as a function of shield thickness at prob- 
abilities a = 0. 001 for a spread of means from approximately 50 percent to 
the 99. 9-percent level, as compared to other researchers at the same a. 
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Figure 9. Solar proton dose curves as a function of shield thickness at prob- 
abilities cl = 0. 01 for a spread of means from approximately 50 percent to the 
99. 9-percent level, as compared to other researchers at the same a. 
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Figure 10. Solar proton dose curves as a function of shield thickness at prob- 
abilities a = 0. 10 for a spread of means from approximately 50 percent to the 
99. 9-percent level, as compared to other researchers at the same a. 
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Figure 11. 99-percent confidence level doses as a function of shield thickness 
for the quiet period of the solar cycle using X = 0. 02 event/week. 
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99. O-percent column of Table 8 to ascertain the number of events for the time 
period depicted. The dose data, D , from model II of Table 3 is used. This 

O 

corresponds to the small events of the 19th cycle, as shown in Table 2. It 
should be recalled that if the reader wishes to estimate the realistic dose to 
an astronaut in a complex spacecraft, then the doses of Figure li should be 
reduced by a factor of about three. This also applies to Figures 7 through 10. 

SECTION V. A NEW PROBABILITY MODEL IS FOUND 


The previous section seemingly provides some rationale for choosing 
a probability function for a given set of conditions: (1) the active 300 weeks 
of a very active cycle such as the 19th, (2) the active 300 weeks of an average 
cycle, and (3) the quiet 270 weeks of any solar cycle. However, if one ex- 
amines Tables 12 through 15, a question arises. If a very high confidence 
level is chosen for the Poisson mean and then the number of events that would 
give a low probability a of getting a worse situation is determined, it seems 
that a very improbable situation is being discussed. That is, if the mean is 
at the 99. 9 percent confidence level, what is the joint probability of seeing 
more than events at a probability no greater than 0. 001? For example, in 
Table 12, one sees that = 25 events when P = 0. 999 and a - 0. 001. At 

first, one might suspect that the chance of both conditions occurring is on the 
order of 10 -6 ; however, care must be taken since one is dealing with cumula- 
tive distribution functions. Thus, the product (l-P)o! does not correspond to 
a unique value of N . This can be ascertained by examining the values of 

(l-P)o! corresponding to = 13 in Table 12. Now, the question is asked: 

"If Xq events are observed in a given time T^, what is the prob- 
ability of seeing more than N events in any observation time t?" 

With this knowledge, a unique value of the probability of seeing N events can 
be made for any period t ^ T Q . 

This question is also important because it is proposed that the values 
of A calculated for the total period of T^ = 300 weeks be used and that this 

value be applied to periods of T . Thus , one might be suspicious of the 
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true probability when one applies the value of A to a time period of, for ex- 
ample, only 50 weeks. However, the reader surely agrees that the best esti- 
mate of X corresponds to the total sample space of the 300 most active weeks 
during the 19th solar cycle. 

To provide an answer to the above, one must first ask, what is the 
distribution of possible Poisson means if in a time there are only 

events observed? If equation (43) , which gives the cumulative distribution of 
P as a function of the upper bound values of m for a given X ^ , is examined, 

one sees that the derivative ( yields the desired upper bound probability 

density function for m , 


-d/3 

dm 


= f(m) 


X„ 

-m 0 
e m 

x o ! 


(46) 


This function is continuous in m and is the so-called gamma distribution in 
the variate m. The expected value of m is E(m) = X + i , the variance is 
a 2 (m) = X Q + 1, and 0 


OO 


/ 


X n 

-m 0 , 
e m dm 


( x o 


V 


‘) 


= i 


(47) 


To summarize the meaning of equation ( 46) , one can state that for a 
given observation of X^ events, the probability of m being in the interval 

m to m + dm is given by the probability equation, 

X„ 

-m 0 , 

. e m dm . 

f(m)dm = — — : . (48) 

V 

See Figure 12 for an illustration of equation (48) with X^ = 3. 

However, it is desirable to investigate the more subtle relationship 
that reflects the distribution of A's for a given observation of X^ events 

over a period of T^ weeks. Therefore, one makes the change of variables 

denoted by m = AT^; dm = T^dA, and our probability density function in the 

variate X becomes 
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m 


Figure 12. Plot of the continuous gamma probability density function when = 3. 


f (A) dA = 


dA 


(49) 


x + i 

T o x o “V 

A e 


V 


For a given value of X , the probability of seeing exactly x events in time t 
is given by the discrete Poisson distribution function, 


, ,x -At 

P(x) = Pr(x; At) = 7~r , x = 0, i, 2, . 


x] 


(50) 


where the probability of having a value of A in the interval A + dA is given by 
the density function of equation (49) and t — T (actually the case t > T is 
equally valid) . 

Using the above definitions, the relationship that is sought is given by 
the following: 


Pr(x,t|X 0 ,T 0 ) = / Pr(x; At) • f(A)dA 


(51) 


Thus, it is being stated that for equation (50) of the Poisson distribution, the 
probability of each possible A (a spectrum of possible means) be folded into 
the equation and the results integrated over all possible values of A from 
zero to infinity. The result is a probability density function that is the 
expected value of the Poisson distribution over all possible means A. Thus, 


X rt + 1 
T„ ° t X 


- . “ X +x -(T +t)A 

Pr M x o- T o)° ,x 0 :-xi- / x e dx - 


(52) 


After integration the following is obtained: 


Pr ( x ’ t l x o’ T o) = 


X n + i 
T 0 ° t X 

V x: 


( X 0 +X ) J 


( T o + t ) 


x o +x + 1 


(53) 
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If the substitution 9 = t/T Q is made and the resulting equation simplified, then 
the discrete distribution function in the variate x is , 


( x + x o) : 
e ) - xfx^r 




x + X 0 + i 


(54) 


(i + 0 ) 

It can be shown that for the discrete distribution function above 


2 ( x+x q) : 


9 


x=0 x!X ! x+X +1 

0 (i + 9 ) 


= 1 


(55) 


and the mean or expected value of x is given by 


x = E(x) = (X Q + i) 9 . (56) 

The proof of equation ( 55) can be shown by using hyper geometric functions. 
The results of equation ( 56) can be found by multiplying the summand of 
equation (55) by x, and after simplifying, the value (X Q + 1)0 can be factored 

out leaving a sum that is the equivalent to that shown in equation (55) . An 
illustration of the density function [equation (54)] is shown in Figure 13 for 
9 = 1 and 9 = 0. 52®. 


It is of some interest to note that when 9 = 1 (t = T Q ) , the function 
becomes simply, 



( x + x o) : 


x I X Q l 2 


x + x o + 1 


(57) 


and 

x = E(x) = X Q + 1 

This report is not intended to be a study of probability theory, and the rami- 
fications of the above distribution function and its possible parallels in other 
statistical work will not be pursued. 


5. 
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The probability distribution of equation (54) is a special form of 
the so called "negative binomial. " 
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Because of the format of Tables 12 through 17 and the manner in which 
probability statements have previously been made , it is desirable to find the 
probability that more than N events are seen, given X Q and 0; or 

y (x + X ) 1 -x 

Pr(x>N|X o> 0) = 1 - Pr(x — N |X Q ,0) = 1 - X=Q x » x » x+X +1 

0 (1 + 0 ) 0 

(58) 


where 0 = — , X = number of events observed in a time T (300 weeks) , and 
*0 ^ 

t is the observation time during which N events are seen. 

Thus, the ultimate relation that is sought to answer the probability 
questions is given by equation ( 58) or the obvious variations associated with it. 

In fact, if the foregoing is valid, one may dispense with the difficulties of 
choosing an upper bound value of the Poisson mean at a given level and then 
determining the probability that N events will be exceeded at some probability 
ot as shown in Tables 12 and 17, where the true probability is actually not known. 

Using equation (58) , the number of combinations of 0, X Q , and N can 

readily become excessive. For this reason, only the values of X^ observed 

events from 0 through 9 will be used with nine different values of 0 from 13 
weeks to 300 (0 = 1) . The N-values will be extended to the point where 
Pr(x>N |X q, 0)<iO 6 . (A summary of these tabulations is shown in Table 20 at 

the end of this section. ) There are many observations that can be made con- 
cerning the comparisons of Table 20 with those of Tables 12 through 15 and 
Table 17. For example, the probability of seeing more than 25 events as 
shown in Table 12 is actually 7. 62 x to” 6 rather than the conjectured 1 x 10 _e 
given at the beginning of this section. Also, the probability of seeing more 
than 10 events when t = 78 weeks with X Q = 3 is 5. 61 x io e , and the probability 

of seeing more than two events is actually 0. 107, which one might guess to be 
less than 0. 100 from Table 17. 

It is very interesting to note that when X Q = 0 and 0 = 1, the chance of 

seeing more than three events is as high as 0. 0625. This infers that even though 
no events were observed during a given 300-week active cycle , one cannot be 
more than 93. 75 percent confident that three or less events could occur in a 
similar cycle. Also, it can be seen from the same column under X Q = 0 that 
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one is 50 percent certain that more than zero events will be seen. This seems 
to infer a dilemma bordering on the naive statement that if you know nothing 
about the probability of an event you can only be 50 percent certain that it 
cannot happen, e. g. , the probability of life on Mars. This last statement 
seems to cast doubt on the usefulness of the case when X Q = 0. 

However, it is of academic interest to investigate further the case when 
X Q = 0. For example, equation (54) becomes 

G X 

Pr(x|0, 6) = ■ ; (59) 

( 1 + 0 ) 

and for 9 > 0, one sees that when x = 0, 

Pr(x = 0 1 0, 6) = -JTJ . (60) 


and 


Pr(x > olo, = l - YTJ = T*~e ' (61) 

The probability of seeing an event as 6 approaches zero becomes very small, 
as one would suspect for a very rare event. As 6 increases to large values 
(0 » 1) , the value of the probability approaches one. This infers that if an 
event can happen at a given small probability, then it is almost certain that the 
event will occur after a sufficiently long period. Thus, equations (60) and (61) 
do not defy intuition in an ordinary sense, but leave us with a rather insecure 
feeling since the number of actual observed events in time T^ is zero. How- 
ever, if one examines the density function of possible values of A [equation 
(49)] when = 0, a plausible probability density function is found, 

-v 

f (A) dX = T e dA . (62) 

Table 18 is presented as a survey of equation (58) for Pr(x > N|3, 6 ) s e 
at several values of e and 0. A point of interest concerns the values of the 
results in Table 17 (Section IV) as compared to those of Table 18. It should 
be clear that the results of Table 18 do not provide the same type of probability 
statements as do those of Table 17. In general, the author feels that the re- 
sults of Table 18 are more useful and to the point; however, others may disagree. 
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TABLE 18. NUMBER OF EVENTS N THAT SATISFY Pr(x>N|3, 9) < e 


t 

(weeks) 

n. e 

tN. 

0. 10 

0. 01 

13 

0. 0433 

1 

2 

26 

0. 0866 

1 

2 

52 

0. 1733 

2 

4 

78 

0. 2600 

3 

5 

104 

0. 3467 

3 

6 

156 

0. 5200 

4 

8 

208 

0. 6933 

6 

9 

260 

0. 8667 

7 

11 

300 

1.000 

9 

13 


0. 001 

3 

4 

5 

7 

8 

10 

13 

15 

17 


0. 0001 

3 

5 

7 

8 

10 

13 

16 

19 

21 


0. 00001 

4 

6 

8 

10 

12 

16 

19 

22 

25 


In any case, one may ascertain which probability level of the parameter X is 
more credible when using the methods of Section IV and constructing tables 
similar to Tables 12 and 17. For example, one sees that the value of X at the 
90-percent level yields event numbers N at probabilities o> in Table 17 that 
seem to be comparable to the probabilities e in Table 18. 


Table 19 is presented for the case when X Q = 0 even though doubt has 

been cast on the validity of the meaning of this rather extreme case. However, 
it does signify a sort of boundary condition for those periods when the prob- 
ability of a large event is very small, as perhaps exists for the active years of 
an average type solar cycle. 


Table 20, which was discussed previously is presented on pages 
58 through 63. 


SECTION VI. REVIEW AND COMMENTS 


In the previous sections , the author has attempted to convey a method 
of considering the radiation hazard associated with solar proton events. The 
method of approach is believed to be more important than the actual results 
presented. The methods and models as developed utilize very little solar 
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TABLE 19. NUMBER OF EVENTS N THAT SATISFY Pr(x>N 1 0, 0)== e 


t 

(weeks) 


0. 1 

0. 01 

0. 001 

0. 0001 

0. 00001 

13 

0. 0433 

0 

1 

2 

2 

3 

26 

0. 0866 

0 

1 

2 

3 

4 

52 

0. 1733 

1 

2 

3 

4 

6 

78 

0. 2600 

i 

2 

4 

5 

7 

104 

0. 3467 

1 

3 

5 

6 

8 

156 

0. 5200 

2 

4 

6 

8 

10 

208 

0. 6933 

2 

5 

7 

10 

12 

260 

0. 8667 

2 

5 

8 

11 

14 

300 

1. 0000 

3 

6 

9 

13 

16 


physics as such and, consequently, will be very unsatisfying to many physicists 
who have examined various aspects of the problem. There is no attempt to 
model or predict a solar proton spectrum, the time dependence of particle 
arrival, or angular distributions of the flux. 

In summary, the total dose in rads-tissue from solar protons was tabu- 
lated for weekly time intervals, and the number of weeks that gave a dose above 
25 rads behind 10 g/cm 2 of aluminum for the active six years of the 19th cycle 
were called dangerous or large event weeks. The number of such event weeks 
was found to be only three weeks during the past 20 years. 

Even though the possibility for smaller events was examined, it was 
found that for any reasonable high confidence level ( 95 percent) , the smaller 
events could be ignored. Consequently, the total particle flux was taken for the 
1 9th cycle , this spectrum was divided by a factor of three , and a single large 
event week was derived. 


The spectrum for this large event is 

J(>p) = 1.76 x 10»e- p/73 (30^ E ^ 115 MeV) 

J(>p) = 3.8x 10 10 e " p / 100 (E >115 MeV) 
where p is in rigidity units of MV. 


(63) 
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TABLE 20. VALUES OF Pr(x > N|X Q , 0) E FOxi A GIVEN 0 =t/T 

0 = i (300 Weeks) 


N \ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

0 

5 . 00 - 1 ” 

7 . 50-1 

8 . 75-1 

9 . 37-1 

9 . 69-1 

9 . 84-1 

9 . 92-1 

9 . 96-1 

9 . 98-1 

1 

2 . 50-1 

5 . 00-1 

6 . 87-1 

8 . 12-1 

8 . 91-1 

9 . 37-1 

9 . 65-1 

9 . 80-1 

9 . 89-1 

2 

1 . 25-1 

3 . 12-1 

5 . 00-1 

6 . 56-1 

7 . 73-1 

8 . 55-1 

9 . 10-1 

9 . 45-1 

9 . 67-1 

3 

6 . 25-2 

1 . 87-1 

3 . 44-1 

5 . 00-1 

6 . 37-1 

7 . 46-1 

8 . 28-1 

8 . 87-1 

9 . 27-1 

4 

3 . 12-2 

1 . 09-1 

2 . 27-1 

3 . 63-1 

5 . 00-1 

6 . 23-1 

7 . 26-1 

8 . 06-1 

8 . 67-1 

5 

1 . 56-2 

6 . 25-2 

1 . 45-1 

2 . 54-1 

3 . 77-1 

5 . 00-1 

6 . 13-1 

7 . 09-1 

7 . 88-1 

6 

7 . 81-3 

3 . 52-2 

8 . 98-2 

1 . 72-1 

2 . 74-1 

3 . 87-1 

5 . 00-1 

6 . 05-1 

6 . 96-1 

7 

3 . 91-3 

1 . 95-2 

5 . 47-2 

1 . 13-1 

1 . 94-1 

2 . 91-1 

3 . 95-1 

5 . 00-1 

5 . 98-1 

8 

1 . 95-3 

1 . 07-2 

3 . 27-2 

7 . 30-2 

1 . 33-1 

2 . 12-1 

3 . 04-1 

4 . 02-1 

5 . 00-1 

9 

9 . 77-4 

5 . 86-3 

1 . 93-2 

4 . 61-2 

8 . 98-2 

1 . 51-1 

2 . 27-1 

3 . 15-1 

4 . 07-1 

10 

4 . 88-4 

3 . 17-3 

1 . 12-2 

2 . 87-2 

5 . 92-2 

1 . 05-1 

1 . 66-1 

2 . 40-1 

3 . 24-1 

11 

2 . 44-4 

1 . 71-3 

6 . 47-3 

1 . 76-2 

3 . 84-2 

7 . 17-2 

1 . 19-1 

1 . 80-1 

2 . 52-1 

12 

1 . 22-4 

9 . 16-4 

3 . 69-3 

1 . 06-2 

2 . 45-2 

4 . 81-2 

8 . 35-2 

1 . 32-1 

1 . 92-1 

13 

6 . 10-5 

4 . 88-4 

2 . 09-3 

6 . 36-3 

1 . 54-2 

3 . 18-2 

5 . 77-2 

9 . 46-2 

1 . 43-1 

14 

3 . 05-5 

2 . 59-4 

1 . 17-3 

3 . 77-3 

9 . 61-3 

2 . 07-2 

3 . 92-2 

6 . 69-2 

1 . 05-1 

15 

1 . 53-5 

1 . 37-4 

6 . 56-4 

2 . 21-3 

5 . 91-3 

1 . 33-2 

2 . 62-2 

4 . 66-2 

7 . 58-2 

16 

7 . 63-6 

7 . 25-5 

3 . 64-4 

1 . 29-3 

3 . 60-3 

8 . 45-3 

1 . 73-2 

3 . 20-2 

5 . 39-2 

17 

3 . 81-6 

3 . 81-5 

2 . 01-4 

7 . 45-4 

2 . 17-3 

5 . 31-3 

1 . 13-2 

2 . 16-2 

3 . 78-2 

18 

1 . 91-6 

2 . 00-5 

1 . 11-4 

4 . 28-4 

1 . 30-3 

3 . 31-3 

7 . 32-3 

1 . 45-2 

2 . 61-2 

19 

9 . 54-7 

1 . 05-5 

6 . 06-5 

2 . 44-4 

7 . 72-4 

2 . 04-3 

4 . 68-3 

9 . 58-3 

1 . 78-2 

20 


5 . 48-6 

3 . 30-5 

1 . 39-4 

4 . 55-4 

1 . 25-3 

2 . 96-3 

6 . 27-3 

1 . 21-2 

21 


2 . 86-6 

1 . 79-5 

7 . 83-5 

2 . 67-4 

7 . 57-4 

1 . 86-3 

4 . 07-3 

8 . 06-3 

22 


1 . 49-6 

9 . 72-6 

4 . 40-5 

1 . 55-4 

4 . 56-4 

1 . 16-3 

2 . 61-3 

5 . 34-3 

23 


7 . 75-7 

5 . 25-6 

2 . 46-5 

9 . 00-5 

2 . 73-4 

7 . 15-4 

1 . 66-3 

3 . 50-3 

24 



2 . 82-6 

1 . 37-5 

5 . 19-5 

1 . 62-4 

4 . 39-4 

1 . 05-3 

2 . 28-3 

25 



1 . 52-6 

7 . 62-6 

2 . 97-5 

9 . 61-5 

2 . 68-4 

6 . 59-4 

1 . 47-3 

26 



8 . 12-7 

4 . 22-6 

1 . 70-5 

5 . 65-5 

1 . 62-4 

4 . 11-4 

9 . 39-4 

27 




2 . 32-6 

9 . 65-6 

3 . 31-5 

9 . 76-5 

2 . 54-4 

5 . 97-4 

28 




1 . 28-6 

5 . 47-6 

1 . 93-5 

5 . 84-5 

1 . 56-4 

3 . 76-4 

29 




7 . 00-7 

3 . 08-6 

1 . 12-5 

3 . 48-5 

9 . 55-5 

2 . 36-4 

30 





1 . 74-6 

6 . 46-6 

2 . 06-5 

5 . 81-5 

1 . 47-4 

31 




1 


3 . 71-6 

1 . 22-5 

3 . 51-5 

9 . 11-5 

32 






2 . 13-6 

7 . 15-6 

2 . 11-5 

5 . 61-5 

33 






1 . 21-6 

4 . 19-6 

1 . 27-5 

3 . 44-5 

34 






6 . 91-7 

2 . 44-6 

7 . 55-6 

2 . 10-5 

35 







1 . 42-6 

4 . 48-6 

1 . 27-5 

36 







8 . 24-7 

2 . 65-6 

7 . 69-6 

37 








1 . 56-6 

4 . 63-6 

38 








9 . 14-7 

2 . 78-6 

39 









1 . 66-6 

40 









9 . 85-7 

41 










42 










43 











9 . 

9 . 

9 . 

9 . 

9 . 

8 . 

7 . 
6 . 
5 . 

5 . 

4 . 

3 . 
2 . 
2 . 
1 . 
1 . 

8 . 

6 . 

4 . 

3 . 

2 . 

1 . 

1 . 

6 . 

4 . 
2 . 
1 . 
1 . 
8 . 

5 . 

3 . 

2 . 

1 . 

8 . 

5 . 

3 . 
2 . 
1 . 
7 . 

4 . 

2 . 

1 . 

1 . 

6 . 


a. Pr(x>N|X o ,0) = i-Pr(x <N|X 0 , 0) 

= Pr(x> N+i|X 0 , 0) 

b. The Notation: 

a. bc-y = a. be x 10~^ 


0 


9 

99-1 

94-1 

81-1 

54-1 

10-1 

49-1 

73-1 

85-1 

93-1 

00-1 

12-1 

32 - 1 
62-1 
02-1 
54-1 
15-1 
43-2 
10-2 
36-2 
07-2 

14 - 2 
47-2 
00-2 
77-3 
52-3 
99-3 
97-3 

28 - 3 

29 - 4 

33 - 4 

40-4 

15 - 4 
36-4 
51-5 

30 - 5 
29-5 
03-5 
25-5 
62-6 
63-6 

81-6 

70-6 

02-6 

16 - 7 


58 


E 


UK -I* 


4 . 64-1 7 . 13-1 8 . 46-1 9 . 18-1 9 . 56-1 9 . 76-1 9 . 87-1 9 . 93-1 9 . 96-1 
2 . 16-1 4 . 47-1 6 . 32-1 7 . 65-1 8 . 53-1 9 . 11-1 9 . 46-1 9 . 68-1 9 . 81-1 


9 4 . 65-4 2 . 96-3 1 . 03-2 2 . 60-2 5 . 35-2 9 . 46-2 1 . 50-1 2 . 17-1 2 . 94-1 


10 2 . 16-4 1 . 49-3 5 . 58-3 1 . 51-2 3 . 29-2 6 . 15-2 1 . 02-1 1 . 56-1 2 . 20-1 2 . 92-1 

11 1 . 00-4 7 . 45-4 2 . 99-3 8 . 61-3 1 . 99-2 3 . 92-2 6 . 86-2 1 . 09-1 1 . 60-1 2 . 22-1 


17 1 . 00-6 1 . 07-5 6 . 00-5 2 . 36-4 7 . 32-4 1 . 90-3 4 . 30-3 8 . 70-3 1 . 61-2 2 . 75-2 

18 4 . 66-7 5 . 22-6 3 . 07-5 1 . 26-4 4 . 07-4 1 . 10-3 2 . 58-3 5 . 42-3 1 . 04-2 1 . 83-2 

19 2 . 54-6 1 . 56-5 6 . 69-5 2 . 25-4 6 . 31-4 1 . 54-3 3 . 34-3 6 . 60-3 1 . 20-2 

20 1 . 23-6 7 . 90-6 3 . 53-5 1 . 23-4 3 . 59-4 9 . 06-4 2 . 04-3 4 . 16-3 7 . 82-3 


9 . 

, 96 - 

-1 

9 . 

, 81 - 

-1 

9 . 

, 46 - 

-1 

8 . 

, 86 - 

-1 

8 . 

, 02 - 

-1 

7 , 

, 01 - 

-1 

5 . 

, 92 - 

-1 

4 . 

, 83 - 

-1 

3 , 

, 82 - 

-1 

2 . 

, 94 - 

-1 

2 , 

, 20 - 

-1 

1 . 

, 60 - 

-1 

1 , 

, 15 - 

-1 

8 , 

, 01 - 

-2 

5 . 

, 49 - 

-2 

3 . 

, 70 - 

-2 

2 . 

, 46 - 

-2 

1 . 

, 61 - 

-2 

1 . 

, 04 - 

-2 

6 i 

, 60 - 

-3 

4 . 

, 16 - 

-3 

2 . 

, 59 - 

-3 

1 , 

, 60 - 

-3 

9 . 

, 74 - 

-4 

5 , 

, 90 - 

-4 

3 , 

, 54 - 

-4 

2 . 

, 11 - 

-4 

1 , 

, 25 - 

-4 

7 , 

, 31 - 

-5 

4 , 

, 26 - 

-5 

2 . 

, 47 - 

-5 

1 , 

, 42 - 

-5 

8 , 

, 14 - 

-6 

4 , 

, 64 - 

-6 

2 , 

, 63 - 

-6 

1 , 

. 49 - 

-6 

8 , 

. 35 - 

-7 


9 . 98-1 

9 . 89-1 

9 . 66-1 

9 . 23-1 

8 . 58-1 

7 . 74-1 

6 . 77-1 

5 . 73-1 

4 . 71-1 

3 . 76-1 


1 . 34-6 

7 . 60-7 











TABLE 20. ( Continued) 

0 = 0.6933 (208 Weeks) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

4 . 09-1 

6 . 51-1 

7 . 94-1 

rH 

CO 

r-s 

« 

CO 

9 . 28-1 

9 . 58-1 

9 . 75-1 

9 . 85-1 

9 . 91-1 

9 . 95-1 

1 . 68-1 

3 . 66-1 

5 . 41-1 

6 . 79-1 

7 . 81-1 

8 . 53-1 

9 . 03-1 

9 . 37-1 

9 . 59-1 

9 . 74-1 

6 . 86-2 

1 . 90-1 

3 . 34-1 

4 . 75-1 

6 . 01-1 

7 . 04-1 

7 . 86-1 

8 . 47-1 

8 . 93-1 

9 . 26-1 

2 . 81-2 

9 . 45-2 

1 . 93-1 

3 . 08-1 

4 . 28-1 

5 . 41-1 

6 . 41-1 

7 . 26-1 

7 . 94-1 

8 . 48-1 

1 . 15-2 

4 . 55-2 

1 . 06-1 

1 . 89-1 

2 . 87-1 

3 . 91-1 

4 . 93-1 

5 . 88-1 

6 . 73-1 

7 . 45-1 

4 . 71-3 

2 . 14-2 

5 . 59-2 

1 . 10-1 

1 . 82-1 

2 . 68-1 

3 . 60-1 

4 . 54-1 

5 . 43-1 

6 . 26-1 

1 . 93-3 

9 . 90-3 

2 . 87-2 

6 . 21-2 

1 . 11-1 

1 . 75-1 

2 . 51-1 

3 . 34-1 

4 . 20-1 

5 . 04-1 

7 . 90-4 

4 . 52-3 

1 . 44-2 

3 . 40-2 

6 . 57-2 

1 . 11-1 

1 . 68-1 

2 . 36-1 

3 . 11-1 

3 . 90-1 

3 . 23-4 

2 . 04-3 

7 . 12-3 

1 . 81-2 

3 . 76-2 

6 . 75-2 

1 . 09-1 

1 . 61-1 

2 . 22-1 

2 . 91-1 

1 . 32-4 

9 . 15-4 

3 . 45-3 

9 . 46-3 

2 . 10-2 

4 . 00-2 

6 . 81-2 

1 . 06-1 

1 . 54-1 

2 . 10-1 

5 . 42-5 

4 . 06-4 

1 . 65-3 

4 . 85-3 

1 . 15-2 

2 . 31-2 

4 . 16-2 

6 . 80-2 

1 . 03-1 

1 . 47-1 

2 . 22-5 

1 . 80-4 

7 . 83-4 

2 . 45-3 

6 . 13-3 

1 . 31-2 

2 . 48-2 

4 . 25-2 

6 . 73-2 

9 . 99-2 

9 . 09-6 

7 . 89-5 

3 . 67-4 

1 . 22-3 

3 . 23-3 

7 . 27-3 

1 . 44-2 

2 . 59-2 

4 . 28-2 

6 . 62-2 

3 . 72-6 

3 . 45-5 

1 . 71-4 

6 . 00-4 

1 . 68-3 

3 . 97-3 

8 . 25-3 

1 . 55-2 

2 . 67-2 

4 . 29-2 

1 . 52-6 

1 . 50-5 

7 . 88-5 

2 . 92-4 

8 . 59-4 

2 . 13-3 

4 . 64-3 

9 . 08-3 

1 . 63-2 

2 . 72-2 

6 . 24-7 

6 . 52-6 

3 . 61-5 

1 . 41-4 

4 . 35-4 

1 . 13-3 

2 . 57-3 

5 . 23-3 

9 . 76-3 

1 . 69-2 


2 . 82-6 

1 . 65-5 

6 . 74-5 

2 . 18-4 

5 . 91-4 

1 . 40-3 

2 . 97-3 

5 . 75-3 

1 . 03-2 


1 . 22-6 

7 . 46-6 

3 . 20-5 

1 . 08-4 

3 . 06-4 

7 . 54-4 

1 . 66-3 

3 . 33-3 

6 . 19-3 


5 . 23-7 

3 . 36-6 

1 . 51-5 

5 . 32-5 

1 . 57-4 

4 . 01-4 

9 . 17-4 

1 . 91-3 

3 . 66-3 



1 . 51-6 

7 . 07-6 

2 . 60-5 

7 . 95-5 

2 . 11-4 

5 . 00-4 

1 . 08-3 

2 . 13-3 



6 . 74-7 

3 . 29-6 

1 . 26-5 

4 . 00-5 

1 . 10-4 

2 . 70-4 

6 . 00-4 

1 . 23-3 




1 . 52-6 

6 . 05-6 

1 . 99-5 

5 . 69-5 

1 . 44-4 

3 . 31-4 

6 . 98-4 




7 . 02-7 

2 . 89-6 

9 . 87-6 

2 . 91-5 

7 . 62-5 

1 . 80-4 

3 . 92-4 





1 . 37-6 

4 . 85-6 

1 . 48-5 

3 . 99-5 

9 . 74-5 

2 . 18-4 





6 . 48-7 

2 . 37-6 

7 . 45-6 

2 . 07-5 

5 . 22-5 

1 . 20-4 






1 . 15-6 

3 . 73-6 

1 . 07-5 

2 . 77-5 

6 . 55-5 






5 . 56-7 

1 . 85-6 

5 . 47-6 

1 . 46-5 

3 . 54-5 







9 . 15-7 

2 . 78-6 

7 . 61-6 

1 . 90-5 








1 . 41-6 

3 . 94-6 

1 . 01-5 








7 . 05-7 

2 . 03-6 

5 . 34-6 


1 . 04-6 2 . 80-6 
5 . 26 - 7 ) 1 . 46-6 
7 . 51-7 







TABLE 20. (Continued) 


Q= 0.5200 (156 Weeks) 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

3 . 42-1 

5 . 67-1 

7 . 15-1 

8 . 13-1 

8 . 77-1 

9 . 19-1 

9 . 47-1 

9 . 65-1 

9 . 77-1 

9 . 85-1 

1 

1 . 17-1 

2 . 71-1 

4 . 23-1 

5 . 56-1 

6 . 66-1 

7 . 52-1 

8 . 19-1 

8 . 69-1 

9 . 06-1 

9 . 33-1 

2 

4 . 00-2 

1 . 19-1 

2 . 23-1 

3 . 37-1 

4 . 50-1 

5 . 53-1 

6 . 44-1 

7 . 21-1 

7 . 84-1 

8 . 35-1 

3 

1 . 37-2 

4 . 97-2 

1 . 09-1 

1 . 87-1 

2 . 77-1 

3 . 71-1 

4 . 65-1 

5 . 52-1 

6 . 32-1 

7 . 01-1 

4 

4 . 69-3 

2 . 01-2 

5 . 05-2 

9 . 72-2 

1 . 59-1 

2 . 31-1 

3 . 11-1 

3 . 94-1 

4 . 75-1 

5 . 52-1 

5 

1 . 60-3 

7 . 93-3 

2 . 25-2 

4 . 81-2 

8 . 59-2 

1 . 36-1 

1 . 96-1 

2 . 63-1 

3 . 36-1 

4 . 10-1 

6 

5 . 48-4 

3 . 07-3 

9 . 72-3 

2 . 28-2 

4 . 44-2 

7 . 56-2 

1 . 17-1 

1 . 67-1 

2 . 25-1 

2 . 88-1 

7 

1 . 88-4 

1 . 18-3 

4 . 10-3 

1 . 05-2 

2 . 21-2 

4 . 04-2 

6 . 65-2 

1 . 01-1 

1 . 43-1 

1 . 93-1 

8 

6 . 42-5 

4 . 44-4 

1 . 69-3 

4 . 71-3 

1 . 07-2 

2 . 08-2 

3 . 65-2 

5 . 85-2 

8 . 75-2 

1 . 24-1 

9 

2 . 20-5 

1 . 66-4 

6 . 89-4 

2 . 06-3 

5 . 01-3 

1 . 04-2 

1 . 93-2 

3 . 27-2 

5 . 15-2 

7 . 61-2 

10 

7 . 51-6 

6 . 19-5 

2 . 76-4 

8 . 88-4 

2 . 30-3 

5 . 08-3 

9 . 96-3 

1 . 77-2 

2 . 93-2 

4 . 53-2 

11 

2 . 57-6 

2 . 29-5 

1 . 10-4 

3 . 76-4 

1 . 03-3 

2 . 42-3 

5 . 00-3 

9 . 36-3 

1 . 62-2 

2 . 61-2 

12 

8 . 79-7 

8 . 40-6 

4 . 30-5 

1 . 57-4 

4 . 57-4 

1 . 13-3 

2 . 45-3 

4 . 81-3 

8 . 70-3 

1 . 47-2 

13 


3 . 07-6 

1 . 67-5 

6 . 47-5 

1 . 99-4 

5 . 16-4 

1 . 18-3 

2 . 42-3 

4 . 57-3 

8 . 02-3 

14 


1 . 12-6 

6 . 46-6 

2 . 64-5 

8 . 54-5 

2 . 33-4 

5 . 56-4 

1 . 19-3 

2 . 35-3 

4 . 29-3 

15 


4 . 05-7 

2 . 48-6 

1 . 07-5 

3 . 62-5 

1 . 03-4 

2 . 58-4 

5 . 79-4 

1 . 18-3 

2 . 25-3 

16 



9 . 43-7 

4 . 26-6 

1 . 52-5 

4 . 54-5 

1 . 18-4 

2 . 76-4 

5 . 87-4 

1 . 15-3 

17 




1 . 69-6 

6 . 31-6 

1 . 97-5 

5 . 34-5 

1 . 29-4 

2 . 86-4 

5 . 83-4 

18 




6 . 66-7 

2 . 60-6 

8 . 44-6 

2 . 38-5 

6 . 00-5 

1 . 37-4 

2 . 90-4 

19 





1 . 06-6 

3 . 59-6 

1 . 05-5 

2 . 74-5 

6 . 50-5 

1 . 42-4 

20 





4 . 27-7 

1 . 51-6 

4 . 58-6 

1 . 24-5 

3 . 04-5 

6 . 85-5 

21 






6 . 28-7 

1 . 98-6 

5 . 54-6 

1 . 40-5 

3 . 27-5 

22 







8 . 47-7 

2 . 45-6 

6 . 42-6 

1 . 54-5 

23 








1 . 07-6 

2 . 91-6 

7 . 18-6 

24 








4 . 65-7 

1 . 30-6 

3 . 31-6 

25 









5 . 76-7 

1 . 51-6 

26 




__ 

L . 


j 



6 . 82-7 


A = 0 . 3467 ( 104 Weeks) 


\x 

N \ 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

2 . 57-1 

4 . 49-1 

5 . 91-1 

6 . 96-1 

7 . 74-1 

8 . 32-1 

8 . 75-1 

9 . 08-1 

9 . 31-1 

9 . 49-1 

1 

6 . 63-2 

1 . 65-1 

2 . 74-1 

3 . 83-1 

4 . 84-1 

5 . 73-1 

6 . 51-1 

7 . 17-1 

7 . 72-1 

8 . 18-1 

2 

1 . 71-2 

5 . 51-2 

1 . 12-1 

1 . 81-1 

2 . 59-1 

3 . 40-1 

4 . 20-1 

4 . 97-1 

5 . 68-1 

6 . 32-1 

3 

4 . 39-3 

1 . 74-2 

4 . 17-2 

7 . 76-2 

1 . 24-1 

1 . 80-1 

2 . 42-1 

3 . 07-1 

3 . 74-1 

4 . 41-1 

4 

1 . 13-3 

5 . 33-3 

1 . 47-2 

3 . 09-2 

5 . 49-2 

8 . 71-2 

1 . 27-1 

1 . 73-1 

2 . 25-1 

2 . 81-1 

5 

2 . 91-4 

1 . 59-3 

4 . 96-3 

1 . 16-2 

2 . 28-2 

3 . 93-2 

6 . 19-2 

9 . 06-2 

1 . 25-1 

1 . 65-1 

6 

7 . 49-5 

4 . 64-4 

1 . 62-3 

4 . 20-3 

8 . 98-3 

1 . 68-2 

2 . 84-2 

4 . 44-2 

6 . 52-2 

9 . 09-2 

7 

1 . 93-5 

1 . 34-4 

5 . 17-4 

1 . 46-3 

3 . 40-3 

6 . 85-3 

1 . 24-2 

2 . 06-2 

3 . 21-2 

4 . 73-2 

8 

4 . 96-6 

3 . 81-5 

1 . 61-4 

4 . 97-4 

1 . 24-3 

2 . 69-3 

5 . 19-3 

9 . 16-3 

1 . 51-2 

2 . 34-2 

9 

1 . 28-6 

1 . 08-5 

4 . 95-5 

1 . 65-4 

4 . 42-4 

1 . 02-3 

2 . 09-3 

3 . 91-3 

6 . 79-3 

1 . 11-2 

10 

3 . 28-7 

3 . 01-6 

1 . 50-5 

5 . 35-5 

1 . 54-4 

3 . 77-4 

8 . 18-4 

1 . 61-3 

2 . 95-3 

5 . 03-3 

11 


8 . 3 8 - 7 ; 

4 . 48-6 

1 . 71-5 

5 . 22-5 

1 . 36-4 

3 . 11-4 

6 . 47-4 

1 . 24-3 

2 . 22-3 

12 



1 . 32-6 

5 . 38-6 

1 . 74-5 

4 . 79-5 

1 . 16-4 

2 . 53-4 

5 . 06-4 

9 . 46-4 

13 



3 . 86-7 

1 . 67-6 

5 . 73-6 

1 . 66-5 

4 . 21-5 

9 . 63-5 

2 . 02-4 

3 . 94-4 

14 




5 . 10-7 

1 . 86-6 

5 . 65-6 

1 . 50-5 

3 . 59-5 

7 . 87-5 

1 . 60-4 

15 





5 . 92-7 

1 . 89-6 

5 . 27-6 

1 . 32-5 

3 . 00-5 

6 . 34-5 

16 






6 . 23-7 

1 . 82-6 

4 . 74-6 

1 . 12-5 

2 . 47-5 

17 







6 . 17-7 

1 . 68-6 

4 . 14-6 

9 . 43-6 

18 








5 . 80-7 

1 . 49-6 

3 . 54-6 

19 









5 . 28-7 

1 . 31-6 

20 

1 









4 . 72-7 
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TABLE 20. (Continued) 


9 = 0.2600 (78 Weeks) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

2 . 06-1 

3 . 70-1 

5 . 00-1 

6 . 03-1 

6 . 85-1 

7 . 50-1 

8 . 02-1 

8 . 43-1 

8 . 75-1 

9 . 01-1 

4 . 26-2 

1 . 10-1 

1 . 91-1 

2 . 76-1 

3 . 60-1 

4 . 41-1 

5 . 15-1 

5 . 83-1 

6 . 43-1 

6 . 96-1 

8 . 79-3 

2 . 97-2 

6 . 29-2 

1 . 07-1 

1 . 59-1 

2 . 17-1 

2 . 79-1 

3 . 41-1 

4 . 04-1 

4 . 64-1 

1 . 81-3 

7 . 57-3 

1 . 90-2 

3 . 71-2 

6 . 23-2 

9 . 43-2 

1 . 32-1 

1 . 75-1 

2 . 23-1 

2 . 72-1 

3 . 74-4 

1 . 86-3 

5 . 39-3 

1 . 19-2 

2 . 23-2 

3 . 72-2 

5 . 68-2 

8 . 13-2 

1 . 10-1 

1 . 44-1 

7 . 72-5 

4 . 45-4 

1 . 47-3 

3 . 63-3 

7 . 49-3 

1 . 36-2 

2 . 25-2 

3 . 47-2 

5 . 03-2 

6 . 96-2 

1 . 59-5 

1 . 04-4 

3 . 85-4 

1 . 05-3 

2 . 38-3 

4 . 70-3 

8 . 38-3 

1 . 38-2 

2 . 13-2 

3 . 13-2 

3 . 29-6 

2 . 42-5 

9 . 87-5 

2 . 96-4 

7 . 26-4 

1 . 55-3 

2 . 96-3 

5 . 19-3 

8 . 52-3 

1 . 32-2 

6 . 78-7 

5 . 52-6 

2 . 47-5 

8 . 07-5 

2 . 14-4 

; 4 . 89-4 

9 . 98-4 

1 . 86-3 

3 . 24-3 

5 . 30-3 


1 . 25-6 

6 . 10-6 

2 . 15-5 

6 . 12-5 

1 . 49-4 

3 . 24-4 

i 6 . 42-4 

1 . 18-3 

2 . 03-3 


2 . 80-7 

1 . 48-6 

5 . 61-6 

1 . 71-5 

4 . 44-5 

1 . 02-4 

1 

2 . 14-4 

4 . 13-4 

7 . 46-4 



3 . 54-7 

1 . 44-6 

4 . 66-6 

1 . 29-5 

3 . 13-5 

6 . 89-5 

1 . 40-4 

2 . 65-4 




3 . 62-7 

1 . 25-6 

3 . 65-6 

9 . 35-6 

2 . 16-5 

4 . 60-5 

9 . 12-5 





3 . 26-7 

1 . 01-6 

2 . 73-6 

6 . 63-6 

1 . 48-5 

3 . 05-5 






2 . 75-7 

7 . 80-7 

1 . 99-6 

4 . 62-6 

9 . 97-6 








5 . 82-7 

1 . 42-6 

3 . 18-6 









4 . 25-7 

9 . 90-7 


9 = 0.1733 (52 Weeks) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

1 . 48-1 

2 . 74-1 

3 . 81-1 

4 . 72-1 

5 . 50-1 

6 . 17-1 

6 . 73-1 

7 . 22-1 

7 . 63-1 

7 . 98-1 

2 . 18-2 

5 . 90-2 

1 . 07-1 

1 . 61-1 

2 . 18-1 

2 . 77-1 

3 . 36-1 

3 . 93-1 

4 . 47-1 

4 . 99-1 

3 . 22-3 

1 . 15-2 

2 . 55-2 

4 . 55-2 

7 . 10-2 

1 . 01-1 

1 . 36-1 

1 . 74-1 

2 . 14-1 

2 . 56-1 

4 . 76-4 

2 . 10-3 

5 . 56-3 

1 . 15-2 

2 . 03-2 

3 . 22-2 

4 . 76-2 

6 . 62-2 

8 . 81-2 

1 . 13-1 

7 . 04-5 

3 . 70-4 

1 . 14-3 

2 . 66-3 

5 . 26-3 

9 . 25-3 

1 . 49-2 

2 . 25-2 

3 . 22-2 

4 . 41-2 

1 . 04-5 

6 . 35-5 

2 . 22-4 

5 . 82-4 

1 . 27-3 

2 . 45-3 

4 . 29-3 

6 . 98-3 

1 . 07-2 

1 . 56-2 

1 . 53-6 

1 . 07-5 

4 . 19-5 

1 . 22-4 

2 . 92-4 

6 . 11-4 

1 . 15-3 

2 . 02-3 

3 . 30-3 

5 . 12-3 

2 . 26-7 

1 . 77-6 

7 . 70-6 

2 . 46-5 

6 . 40-5 

1 . 45-4 

2 . 94-4 

' 5 . 48-4 

9 . 55-4 

: 1 . 57-3 


2 . 89-7 

1 . 38-6 

4 . 81-6 

1 . 36-5 

3 . 29-5 

7 . 15-5 

11 . 42-4 

12 . 62-4 

I 4 . 55-4 



2 . 44-7 

9 . 16-7 

2 . 78-6 

7 . 24-6 

1 . 67-5 

3 . 52-5 

6 . 87-5 

1 . 26-4 





5 . 55-7 

1 . 54-6 

3 . 79-6 

8 . 43-6 

1 . 73-5 

3 . 34-5 






3 . 19-7 

8 . 30-7 

1 . 95-6 

4 . 23-6 

8 . 53-6 








4 . 38-7 

9 . 98-7 

2 . 11-6 

5 . 05-7 
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TABLE 20. (Concluded) 


e= 0.0867 (26 Weeks) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

7 . 98-2 

1 . 53-1 

2 . 21-1 

2 . 83-1 

3 . 40-1 

3 . 93-1 

4 . 41-1 

4 . 86-1 

5 . 27-1 

5 . 64-1 

6 . 36-3 

1 . 81-2 

3 . 42-2 

5 . 41-2 

7 . 69-2 

1 . 02-1 

1 . 29-1 

1 . 58-1 

1 . 87-1 

2 . 17-1 

5 . 07-4 

1 . 91-3 

4 . 49-3 

8 . 44-3 

1 . 39-2 

2 . 09-2 

2 . 96-2 

3 . 98-2 

5 . 15-2 

6 . 47-2 

4 . 05-5 

1 . 89-4 

5 . 32-4 

1 . 16-3 

2 . 18-3 

3 . 67-3 

5 . 74-3 

8 . 45-3 

1 . 19-2 

1 . 61-2 

3 . 23-6 

1 . 81-5 

5 . 91-5 

1 . 47-4 

3 . 09-4 

5 . 77-4 

9 . 89-4 

1 . 58-3 

2 . 41-3 

3 . 50-3 

2 . 57-7 

11 . 68-6 

6 . 25-6 

1 . 75-5 

4 . 07-5 

8 . 35-5 

1 . 56-4 

2 . 70-4 

4 . 40-4 

6 . 84-4 


1 . 52-7 

6 . 37-7 

1 . 98-6 

5 . 07-6 

1 . 13-5 

2 . 28-5 

4 . 25-5 

7 . 42-5 

1 . 23-4 




2 . 13-7 

6 . 01-7 

1 . 45-6 

3 . 16-6 

6 . 30-6 

1 . 17-5 

2 . 06-5 






1 . 76-7 

4 . 14-7 

8 . 82-7 

1 . 75-6 

3 . 25-6 









2 . 44-7 

4 . 82-7 


9 - 0.0433 (13 Weeks) 



0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

4 . 15-2 

8 . 13-2 

1 . 19-1 

1 . 56-1 

1 . 91-1 

2 . 25-1 

2 . 57-1 

2 . 88-1 

3 . 17-1 

3 . 46-1 

1 

1 . 73-3 

5 . 03-3 

9 . 79-3 

1 . 59-2 

2 . 31-2 

3 . 15-2 

4 . 09-2 

5 . 11-2 

6 . 22-2 

7 . 40-2 

2 

7 . 16-5 

2 . 78-4 

6 . 73-4 

1 . 30-3 

2 . 21-3 

3 . 43-3 

4 . 98-3 

6 . 90-3 

9 . 20-3 

1 . 19-2 

3 

2 . 98-6 

1 . 44-5 

4 . 17-5 

9 . 41-5 

1 . 82-4 

3 . 17-4 

5 . 11-4 

7 . 76-4 

1 . 13-3 

1 . 57-3 

4 

1 . 24-7 

7 . 17-7 

2 . 42-6 

6 . 23-6 

1 . 35-5 

2 . 61-5 

4 . 62-5 

7 . 66-5 

1 . 20-4 

1 . 80-4 

5 



1 . 34-7 

3 . 88-7 

9 . 34-7 

1 . 98-6 

3 . 82-6 

6 . 84-6 

1 . 15-5 

1 . 86-5 

6 






1 . 42-7 

2 . 94-7 

5 . 66-7 

1 . 02-6 

1 . 75-6 

7 










1 . 56-7 
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Using this spectrum, one can calculate the tissue dose in rads at the 
center of an aluminum spherical shell (see Table 3) . To correct for geometric 
effects and self- shielding, this dose should be reduced by a factor of approxi- 
mately three. If the space mission is planned during the quiet period of a 
cycle, then the small event dose curve (D3) of model II (Table 3) may be used 
with X = 0. 02 (m = 1 event/year). 

To predict the probability of an event occurring, the Poisson distribu- 
tion seemed to be the most logical choice. Section m was devoted to examin- 
ing this conclusion and the methods of using the probability model. Section IV 
defined the confidence one could use in employing the Poisson process and 
specifically arriving at confidence levels for the experimental or observed 
value of the mean of the Poisson distribution function. Several examples were 
given for different mission lengths, and comparisons were made to the results 
of other authors. 


Section V was an extension of the Poisson process to incorporate the 
concept of small sample theory and arrive at the expected distribution function 
that answers the following question, If events are observed in time T^, 

what is the probability of seeing X events in any observation time t? The re- 
sults are represented by the discrete probability density function in the variate 
x, 


X + i 

T , ° 

X f x T 

Aq . x . 


( T o + t ) 


Table 20, which is presented on pages 58 through 63, was generated using the 
preceding function. 


In the beginning of this work, the author intended to avoid commenting 
on the methods used by other investigators, but in order to explain the radical 
differences shown in Figures 8, 9, and 10, the following comments are per- 
tinent. One of the more common procedures used is to obtain the logarithms 
to the base ten of the solar proton flux above 30 MeV for each of the events of 
the 19th cycle. These data are then plotted on normal probability graph paper, 
and a distribution called the log normal is obtained. 

Thus, if x = log^<£, one has the log normal probability density function, 
, - (x-x) l /2 cr 2 

f(x)dx = — e x dx , (65) 

n/^°x 
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where x and <r 2 are the mean and variance respectively and are estimated 

X 

from the sample variables xi = log^^i. This is a true normal distribution 

in the variate x and there is no virtue in examining the transformed distribu- 
tion in the variable <p, which has a more complex representation. Also, the 
values of the mean and standard deviation between the two distributions are not 
simply antilog related. For example, 

cr ^ 2 

log 10 ? = l°g i0 ? + 2 Io g i0 e (66a) 

and 

°x ! " *°%0 elog 10 ( 1+ -J5-) • < 66W 

A 

where cr^ is the variance in the distribution <p. For more details, see 

Reference 11. The author has assumed that the users of this distribution keep 
their statistics in the variable x, which is certainly the simplest process. 


From Reference 10, 60 events were taken that had fluxes measured 
above 30 MeV. The 60 events chosen are shown in Table 21. The logarithms 
of these 60 entries were ordered from the smallest to the largest, and the 
normalized cumulative sums were plotted on log normal paper as shown in 
Figure 14. The ungrouped data had a mean of 7. 39, and the standard deviation 
was 0. 97 as shown at bottom of Table 21. The straight line in Figure 14 de- 
picts the cumulative normal distribution with a mean of 7. 39 and a variance 
of 0. 97. From an examination of this fit, one might conclude that the log 
normal gives a reasonable representation of the data. Note that the flux for 
the November 12, 1960, event is plotted as a square that falls exactly on the 
line of best fit. This occurs because when the cumulative distribution from 
a discrete set of data is obtained, the last point has the cumulative probability 
of 1. 00; but since this cannot theoretically occur, the best choice is to place 
this point on the best fit line. This point, which is log^^ = 9. 9562, corresponds 
to the cumulative probability of 99. 6 percent, or in terms of the standard 
deviation, the November 12, 1960, event is at 2. 63 standard deviations above 
the mean based on the values shown below Table 21. If one wished to go to the 
99. 9-percent level, it would be necessary to take 3. 09 standard deviations 
above the mean, or a log flux of 10. 4016. If the antilog of this value is found, 
one sees that the flux at the 99. 9-percent level in the log normal distribution 
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TABLE 21 . SOLAR FLARE EVENTS OF THE 19TH CYCLE 


Event 

Flux 

Date 

(E>30 MeV) 

2-23-56 

1 . 00(+9) 

3-10-56 

1.10(+8K 

8-31-56 

2. 50(+7) 

11-13-56 

1. 00(+8) 

1-20-57 

3.00(+8H 

4- 3-57 

5. 00(+7) 

4- 6-57 

3. 80(+7) 

6-21-57 

1. 50(+8) M 

7- 3-57 

2. 00(+7) 

7-24-57 

7. 50(+6) 

8- 9-57 

1. 50(+6) 

8-29-57 

1.20(+8)-/ 

8-31-57 

8. 00(+7) 

9- 2-57 

5. 00(+7) 

9-12-57 

6. 00(+6) 

9-21-57 

1. 50(+6) 

10-20-57 

5. 00(+7) 

11- 4-57 

9. 00(+6) 

2- 9-58 

5. 00(+6) 

3-23-58 

2. 50(+8) t 

3-25-58 

6. 00(+8)t 

4-10-58 

5. 00(+6) 

7- 7-58 

2. 50(+8) , 

7-29-58 

8. 50(+6) 

8-16-58 

4. 00(+7) 

8-22-58 

7. 00(+7) 

8-26-58 

1. 10(+8) : 

9-22-58 

6. 00(+6) 

2-13-59 

2. 80(+7) 

5-10-59 

9. 60(+8) ^ 


Summary of Statistics: 


Logio <t> 

Event 

(E>30 MeV) 

Date 

9.0000 

6-13-59 

8.0414 

7-10-59 

7.3979 

7-14-59 

8.0000 

7-16-59 

8.4771 

8-18-59 

7.6990 

9- 2-59 

7.5798 

1-11-60 

8.1761 

3-29-60 

7.3010 

3-30-60 

6.8751 

4- 1-60 

6.1761 

4- 5-60 

8.0792 

4-28-60 

7.9031 

4-29-60 

7.6990 

5- 4-60 

6.7782 

5- 6-60 

6.1761 

5-13-60 

7.6990 

6- 1-60 

6.9542 

8-12-60 

6.6990 

9- 3-60 

8.3979 

9-26-60 

8.7782 

11-12-60 

6.6990 

11-15-60 

8.3979 

11-20-60 

6. 9294 

7-11-61 

7.6021 

7-12-61 

7.8451 

7-15-61 

8.0414 

7-18-61 

6.7782 

7-20-61 

7. 4472 

7-28-61 

8. 9823 

9- 8-61 


60 

Y x. = 443.5541 

iti 1 

60 

T x. 2 = 3335.90112251 

1=1 1 


x = 7.3926 

a 2 = 0.94828 ; cr 

x x 


Flux 

(E>30 MeV) 

8. 50(+7) 

1. 00(+9) 

1 . 30(+9) 

9. 10(+8)y 
1 . 80(+6) 

1 . 20(+7) 

4. 00(+5) 

6 . 00 (+ 6 ) 
6 . 00 (+ 6 ) 

5. 00(+6) 

1 . 10(+6) 

5. 00(+6) 

7. 00(+6) 

6 . 00 (+ 6 ) 

4. 00( + 6) 

4. 00(+6) 

4. 00(+5) 

6. 00(+5) 

3. 50(+7) 

2 . 00 (+ 6 ) 

9. 04(+9) 

7. 20(+8) 

4. 50(+7) 

3. 00(+6) 

4. 00(+7) 

1. 30(+7) 

3. 00(+8) 

5. 00(+6) 

4. 40(+6) 

3. 00(+6) 


0. 9738 


Log 10 <p 
(E>30 MeV) 

7. 9294 
9.0000 
9.1139 
8.9590 
6.2553 

7. 0792 

5.6021 

6.7782 

6. 7782 

6.6990 

6.0414 

6.6990 

6.8451 

6.7782 

6.6021 

6.6021 

5.6021 

5.7782 
7.5441 
6.3010 

9.9562 

8.8573 

7.6532 

6.4771 

7.6021 

7. 1139 

8.4771 

6.6990 
6.6435 

6.4771 
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100% x PROBABILITY OF FLUX LESS THAN X 
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X= LOG 10 4> 


Figure 14. Cumulative distribution of solar proton event log fluxes 

during solar maximum. 



is 2. 52 x lo 10 ^ > ^ > which is a factor of about 2. 8 larger than 

the November 12, i960, event. Hence, it would seem that a flux above 30 
MeV, which is three times larger than the November 12, 1960, event, would 
have less than a 0. 001 probability of occurring. For example, the probability 
of an event that is ten times larger than the November 12, 1960, event is about 
10~ 4 , or 99. 99 percent of all events would be less than such an event. 

From the above analysis, one may be convinced that a reasonable upper 
bound value for a single event is at most a factor of three larger than the 
November 12, 1960, event. For this reason, the author feels that the results 
obtained by some investigators for the extreme probability tails must depict 
a smaller probability than the estimates given. These investigators have 
possibly used joint probabilities of flux and rigidity parameters, which may be 
a factor of ten or more smaller than those indicated by the 0. 1-percent prob- 
ability tail shown in the various reports at the disposal of the author. For 
example, from Table 18 of this report, one observes that the actual probability 
of seeing more than seven large events during 78 active weeks of the 19th 
cycle is less than 10 -3 . However, in Figure 8, the work of Snyder at 20 g/cm 2 
indicates a dose that requires 48 large events [as defined by equation (63)] at 
the 0. 001 probability level. This seems high unless Mr. Snyder uses a large 
event that is about seven times larger than that used by the author. There is 
no intention here of isolating Mr. Snyder's work [9], since the values of Mr. 
Yucker of McDonnell Douglas [10] are obviously higher than any other pre- 
vious analysis and would require 120 large events, or else each event is 17 
times larger than the largest event of the 19th cycle if only 7 large events are 
seen according to Table 18. (If the reader wishes to use the 99. 9-percent 
probability event in the present work, he may ignore self-shielding and geo- 
metry factors, a factor of three. ) 

Of course, all cumulative distributions have the same general appear- 
ance and tend to mask any special structure of the basic frequency distribu- 
tion. The author, being somewhat perplexed by this state of affairs, has 
attempted to examine the data of Table 21 from a somewhat more fundamental 
approach. These data are grouped into intervals and frequency tables as 
shown in Table 22. The mean and the standard deviation of this grouped data 
are not exactly the same as in the ungrouped data, and for simplicity of 
plotting, the frequency histograms are transformed to units of t = (x- x)/cr. 

The value of a = 1. 0 is used for convenience. This would lead to an error of 
from 1 percent to approximately 4 percent in the correct value of the abscissa, 
but this would not be detected in the plots shown in Figures 15 through 17. The 
ordinates or relative frequencies are shown in the fifth column of the tables 
composing Table 22. The smooth curve drawn on Figures 15 through 17 is 
that of the Gaussian density function, 
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TABLE 22. ANALYSIS OF THE DATA IN TABLE 21 


Class 

Interval, 

Mid 

Frequency 

f. 

i 


f. 

i 


Cumulative Relative 

Log 0 


X 

t=x-7. 4 

At x 60 

End t 

Frequency 

5.6 - 

6.0 

5.8 

3 

-1.6 

0. 12500 

-1.4 

0. 0500 

6.0 - 

6.4 

6.2 

5 

-1.2 

0.20825 

-1.0 

0.1333 

6.4 - 

6.8 

6.6 

15 

-0.8 

0.62500 

-0.6 

0.3833 

6.8 - 

7.2 

7.0 

6 

-0. 4 

0.2500 

-0.2 

0.4833 

7.2 - 

7.6 

7.4 

5 

0 

0.20825 

+ 0.2 

0. 5666 

7.6 - 

8.0 

7.8 

10 

+0.4 

0. 41667 

+ 0.6 

0.7333 

8.0 - 

8.4 

8.2 

6 

+ 0. 8 

0.2500 

+ 1.0 

0.8333 

8.4 - 

8.8 

8.6 

3 

+ 1.2 

0.1250 

+ 1.4 

0.8833 

00 

0O 

1 

9.2 

9.0 

6 

+ 1.6 

0.2500 

+ 1.8 

0. 9833 

9.2 — 

9.6 

9.4 

0 

+ 2.0 

0 

+2.2 

0. 9833 

9.6 - 10.0 

9.8 

1 

+2.4 

0.041667 

+2.6 

1.0000 



2 

60 

1 




i 

Summary of 

Statistics: 









X 

= 7.39, cr — 0. 97 




Class 
Interval, 
Log 0 

Mid 

x 

Frequency 

f. 

i 

t=x-7. 4 

f. 

i 

At x 60 

End t 

Cumulative Relative 
Frequency 

5.6 — 6.4 

6.0 

8 

-1. 4 

0. 1667 

-1.0 

0. 1333 

6.4- 7.2 

6.8 

21 

-0.6 

0.4375 

-0.2 

0. 4833 

7.2 — 8.0 

7.6 

15 

+0.2 

0.3125 

+0.6 

0.7333 

8.0 — 8.8 

8.4 

9 

+ 1.0 

0.1875 

+ 1. 4 

0.8833 

8.8— 9.6 

9.2 

6 

+ 1.8 

0.1250 

+2.2 

0.9833 

9.6-10.4 i 

10.0 

1 

+2.6 

0.0208 

+3. 0 

1.0000 


2 

60 





Summary of Statistics: 








x - 7.42, cr - 0.99 




Class 

Interval, 

Mid 

Frequency 

f. 

l 


f 

X 


Cumulative Relative 

Log 0 

X 

t=x-7. 5 

At x 60 

End t 

Frequency 

5.6 — 6.6 

6.1 

10 

-1. 4 

0. 1666 

-0.9 

0.1667 

6.6— 7.6 

7.1 

24 

-0.4 

0.4000 

+0.1 

0.5667 

7.6 — 8.6 

8.1 

18 

+ 0.6 

0.3000 

+ 1.1 

0. 8667 

8.6 - 9.6 

9.1 

7 

+ 1.6 

0. 1167 

+2.1 

0. 9834 

9.6 - 10.6 

10. 1 

1 

+2.6 

0.0167 

+3.1 

1.0000 


2 

60 





Summary of 

Statistics: 








x = 7.51, a - 0.96 
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Figure 15. Normal density function and grouped frequency data 
from top array of Table 22. 




t 


Figure 16. Normal density function and grouped frequency data 
from middle array of Table 22. 
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Figure 17. Normal density : 
from bottom 


, NORMAL DENSITY FUNCTION 



n and grouped frequency data 
of Table 22. 


1 


(67) 


f (t)dt = 


\l 2ir 



The midpoint of the frequency histograms is plotted at the value of t shown in 
Table 22. The width of these histograms is At, or the simple difference in the 
value of the tabulated t’s. By choosing a value of a = i. 0, At becomes 
equal to the grouping intervals in the variate log^<£ as shown in Table 22. 

This small inaccuracy allows a simple direct comparison between the three 
curves shown in Figures 15 through 17 and the true Gaussian frequency dis- 
tributions. The results shown in Figures 15 through 17 indicate a reasonable 
resemblance to a Gaussian distribution. However, from observation, one 
sees that the distribution has a skew to the right. A measure of skewness 
could be calculated from the data in Tables 21 and 22 but would be of no more 
value than the above observation. The purpose of the preceding is not to find 
the best probability density function to fit the data but to demonstrate that 
the use of the log normal for the particle flux above 30 MeV is reasonable. 
However, a reasonable use of the above log normal distribution would lead 
one to choose a maximum flux (above 30 MeV) that does not exceed the largest 
observed event by more than a factor of three to arrive at the 99. 9-percentile 
level. 


From Table 22 (first data group) , Figure 18 is presented to compare 
the cumulative distribution of the grouped data to the true Gaussian drawn as 
a solid line. It should be noted how well the data follow the Gaussian dis- 
tribution where the cumulative area of the rectangles in Figure 15 is plotted 
at the end points , as they should be to obtain the comparable values on the 
Gaussian ogive. The reader can recognize the difference in implications 
between Figure 18 and Figure 15. Stated simply, a cumulative distribution 
hides the true nature of the goodness of fit of basic data to a given distribution 
function. 

In conclusion, one graph will be presented for the dose expected on 
the skin of an astronaut in a typical spacecraft for a 78-week mission during 
the active weeks of an average cycle at various probabilities, €, as shown 
in Table 19 of Section V. The results are shown in Figure 19. The doses 
in Figure 19 have been made a factor of three lower than a point dose at the 
center of a sphere. Figure 19 presents what the author believes to be a 
realistic evaluation of the risk of solar proton events to space travel. Thus , 
if one assumes that 6. 5 g/cm 2 of aluminum is equal to 5 cm of tissue and the 
spacecraft shield is 13. 5 g/cm 2 , then the 5-cm depth dose is about 17 rads at 
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100% x PROBABILITY OF FLUX LESS THAN X 



Figure 18. Cumulative distribution of grouped data 
from top array of Table 22. 



RADS-TISSUE TO ASTRONAUT SKIN 



Figure 19. Expected dose at probabilities e for a 78-week mission 
on astronaut skin inside typical spacecraft of uniform wall thick- 
ness, x g/cm 2 , during the active part of an average cycle. 
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the 99-percent confidence level during the active part of an average cycle. If 
one believes that the galactic cosmic ray (GCR) dose is about 12 rads/year 
and is unattenuated by the shield (this is in doubt for the high Z component)* 
one sees that in 18 months the OCR dose is 18 rads. Using Hie above assump- 
tions, the GCR depth dose is about the same as the solar proton depth dose for 
an 18-month trip at a 99-percent confidence level during an average solar 
cycle (not the 19th cycle). 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, 35812 
November 30, 1970 124-09-21-00-62 
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